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ABSTRACT 

We present results on low-resolution mid-infrared (MIR) spectra of 70 infrared-luminous galaxies obtained with the 
Infrared Spectrograph (IRS) onboard Spitzer. We selected sources from the European Large Area Infrared Survey 
(ELAIS) with 5i5 > 0.8 mjy and photometric or spectroscopic z > 1. About half of the sample are QSOs in the 
optical, while the remaining sources are galaxies, comprising both obscured AGN and starbursts. Redshifts were 
obtained from optical spectroscopy, photometric redshifts and from the IRS spectra. The later turn out to be reliable 
' for obscured and/or star-forming sources, thus becoming an ideal complement to optical spectroscopy for redshift 

5H ' estimation. 

We estimate monochi'omatic luminosities at several restframe wavelengths, equivalent widths and luminosities 
for the PAH features, and strength of the silicate feature in individual spectra. We also estimate integrated 8-1000 
|im infrared luminosities via spectral energy distribution fitting to MIR and far-infrared (FIR) photometry from the 
Spitzer Wide-area Infrared Extragalactic survey (SWIRE) and the MIR spectrum. Based on these measurements, we 
classify the spectra using well-known infrared diagnostics, as well as a new one that we propose, into three types of 
source: those dominated by an unobscured AGN, mostly corresponding to optical quasars (QSOs), those dominated 
by an obscured AGN, and starburst-dominated sources. Starbursts concentrate at z ~ 0.6-1.0 favored by the shift of 
the 7.7-|xm PAH band into the selection 15-|im band, while AGN spread over the 0.5 < z < 3.1 range. 

Star formation rates (SFR) are estimated for individual sources from the luminosity of the PAH features. An 
estimate of the average PAH luminosity in QSOs and obscured AGN is obtained from the composite spectrum of all 
sources with reliable redshifts. The estimated mean SFR in the QSOs is 50-100 Moyr"' , but the implied FIR lumi- 
nosity is 3-10 times lower than that obtained from stacking analysis of the FIR photometry, suggesting destruction of 
the PAH caniers by energetic photons from the AGN. The SFR estimated in obscured AGN is 2-3 times higher than 
in QSOs of similar MIR luminosity. This discrepancy might not be due to luminosity effects or selection bias alone, 
but could instead indicate a connection between obscuration and star formation. However, the observed correlation 
between silicate absorption and the slope of the near- to mid-infrared spectrum is compatible with the obscuration of 
the AGN emission in these sources being produced in a dust torus. 
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1 INTRODUCTION 

Infrared-luminous Galaxies (LIRGs), those objects with infrared 
(IR, [8-lOOO^m]) luminosities above lO" L© (and the Ultralumi- 
nous Infrared Galaxies, those objects with Lir > 10'- Lq) are an 
important population for understanding the cosmological evolution 
of galaxies. At low redshifts LIRGs are relatively rar e, with only 
a few hundred known over the whole sky at z < 0.2 l lSoifer et al.l 
Il984h . At higher redshifts they are much more common; at z > 1 
they reach a density on the sky of several hundred per square de- 
gree, and make a substantial, perhaps even dominant contribution 
to both the Cosmic Infrare d Background and the volume-averaged 
cosm ic star for mation rate dLagache et al .l2005l : ISanders & Mirabell 
[l996 : Lonsdal eiraLlbOod) . 

Theoretical and empirical studies have shown that many 



acting svstems jGallimore & Keell 19931: C 


ements et alj 19961: 


ISanders& Mirabell 1 19961; Borne et alj 19991: 


Farrahetal.1 12001). 


The power source behind the IR emission was initially controver- 



sial, with both star formation and AGN activity proposed as the 
sole driver, but in recent years a consensus has started to emerge, 
with most studies finding that both starbursts and AGN power 
low-redshift LIRGs and ULIR Gs, with the star for mation usually 
domi nating (Genzel et al. 199 i iFarrah et al]|2003l ; [imanishi et al., 
l2007h . At high redshift s a broadly similar picture is thought to 



exist (Iparrah et al."2002'; 'Chapman et al.ll2003l; iTakata et aljhood : 
iBerta et al. 2007; Bridge et al. 2007). 



The ability to reliably distinguish between star formation and 
AGN activity in high redshift LIRGs is essential for tracing both 
the cosmic star formation history and the evolution of AGN ac- 
tivity as a function of redshift. Since these galaxies are dusty, the 
usual optical/near-infrared diagnostics are not always reliable. In- 
stead, mid-infrared (MIR) spectroscopy has proved itself the key 
tool for determining the dominant energy so urce in dusty galax- 
ies (e.g. [O enzel et al. 1998; Rigopoulou et al. 19991 ; [Laurent et al] 
l200d ; riVan et aU2001l ;IWeedman et al. 2005), but until recently the 
limited sensitivity of the Infrared Space Observatory {ISO) con- 
strained the spectroscopic study of U LIRGs to the local Universe. 
With the Infrared Spectrograph (IRS. lHouck et alj|2004l) onboard 
Spitzer, wavelength coverage and sensitivity have improved enough 
to allow redshift determination and s pectroscopic classification for 
ULIRGs at all redshifts up to z ~ 3 faouck et alj|2005l: 1^ et all 



l2005l ; IWeedman et al.l2006l ; lArmus et alj2007l ; lFarrah et al.l20o'7ir 

In most cases, spectroscopic surveys of high-redshift sources 
with the IR S have selected sou r ces with a high infrared-to-optical 
flux ratio faouck et al.l |2005| ; lYan et alj l2005l : IWeedman et"ai] 
l2006l) . These turn out to be mostly obscured AGN at high redshift, 
with spectra dominated by silicate absorption. Starburst ULIRGs 
with strong PAH emission are rare in these samples, but prolific if 
the selection criteria instead include sub-mm detections iLutz et al.. 



2005 ) or prominent stellar 'bumps' in the IRAC bands jparrah et al.l 
20081) . 

In contrast, the IRS has undertaken relatively few surveys of 
sources selected solely on the basis of MIR flux. This means that 
moderately to lightly absorbed sources will be under represented in 
IRS surveys, thus drawing an incomplete picture of the nature of the 
infrared galaxy population at high redshifts. In this work we present 
results on a sample of 70 LIRGs and ULIRGs selected in the MIR 
with no constraints on the optical fluxes (other than the requirement 
of detection in at least 3 optical bands in deep imaging and photo- 
metric or spectroscopic z > 1. Our objectives are: (a) to classify the 
sources into starburst-dominated (SB), obscured AGN-dominated 
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Figure 1. Ratio of tlie flux to the r-band flux as a function of r mag- 

nitude for the sources in the ELAIS-IRS sample. Asterisks indicate sources 
classified as QSOs from their optical spectrum or SED, while diamonds 
represent sources classified as galaxies. Triangles represent sources with 
uncertain or unknown optical class due to lack of enough optical data. 



(AGN2) and unobscured AGN-dominated (AGNl), based on their 
MIR spectral properties, and compare this classification with that 
obtained from optical data, (b) To estimate the IR luminosity of the 
sources and the relative contribution of AGN and star formation, 
(c) To obtain star formation rates for the obscured and unobscured 
AGN. 

The paper is structured as follows: §2 describes the sample 
selection; §3 discusses the observations and data reduction; §4 
presents the measurements and diagnostics; §5 compares the re- 
sults to other high-z MIR-selected samples; and §6 summarizes our 
conclusions. We adopt the following cosmology: f2,„ = 0.27, Q./^ = 
0.73, Ho = 71 kms"' Mpc"'. 



2 SAMPLE SELECTION 

The sources were selected from the European La rge Area ISO Sur- 
vey (EL AIS ; loiTver et al JI2O OO : Rowan-Robinso n et al.ll2004h final 
band-merged catalog jRowan-Robinsoni2003t) . The selection cri- 
teria required detection in the ISOCAM LW3 band jVaccari et al.l 
[2005) with 5i5 > 0.8 miy, as well as detection in at least 3 optical 
bands and spectroscopic or photometric redshift z ^ 1. 

The 70 objects that mat ched were observed with the IRS 
jHeman-Caballero et alj|2003) . resulting in one of the largest z <; 
1 samples of MIR spectra of infrared-luminous galaxies to date. 
Table [T] indicates the name and coordinates of these sources. Only 
the first 15 entries of each table are shown here. The complete ver- 
sion of the tables is available as on-line material. 

In addition to the MIR spectra, there is also optical and IR 
photometry available for all sources, as well as optical spectroscopy 
for roughly half of them. 

Por sources in the northern hemisphere (ELAIS Nl and N2 
fields) optical U, g, r, i, Z photometry comes from the Isaac 
Newton Telescope Wide Pield Survey (WPS), with 5<r limiting 
magnitudes 23.4, 24.9, 24.0, 23.2 and 21.9 (Vega) respectively 
dOonzalez-Solares et alj200^ . In the southern hemisphere (ELAIS 
SI) optical photom etry was obtained down to S, V ~ 25, R ~ 24.5 
teerta et al.i2006l) . 

Eight sources (six of them in ELAIS SI) have near-infrared 
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(NIR) photometry fro m the Two Micron All Sky Survey (2MASS; 
ISkrutskie et al.ll200^ . and ten more are detected by the UK IRT In- 
frared Deep Sky Survey (UKIDSS: ILawrence et aP 12007 1 in the 
central region of ELAIS Nl. Table [2] summarizes the optical and 
NIR photometry for the ELAIS-IRS sources. 

The three ELAIS fields in which the ELAIS-IRS sample was 
selected were observed by Spitzer in seven IR bands from 3.6 to 
160 |,im, as part of the Spitzer Wide-A rea Infrared Extragalactic 
Survey (SWIRE: .Lonsdale e t Zl l2003al . [2004). We performed a 
source match between the ELAIS-IRS catalog and the latest work- 
ing catalog from SWIRE and found that 69 out of 70 ELAIS-IRS 
sources are detected in at least two IRAC bands and in the 24-|xm 
band from MIPS, while 33 are detected in the 70-|xm band and only 
14 at 160 yim (see TableO. 

One source, EIRS-65, has no infrared counterpart in the 
SWIRE catalog and is undetected in the IRS spectrum. Since the 
I5-[.im detection is strong, the I5-[.im source is probably real but 
was associated to the wrong optical counterpart. 

Photometric redshifts were calculated using optical and NIR 

S ihotometry for all sources in the ELAIS bandmerge catalog 
Rowan-Robinson et alj|2004l) . These estimates were revised with 
posteriori ty to the sample selection usi n g IRA C photometry from 
SWIRE l iRowan-Robinson et alj |2004 l2008h . leading in some 
cases to revised photometric z < 1. At least four and up to nine 
optical and NIR bands are used for photometric redshift estima- 
tion, except in four sources whose redshifts are obtained from just 
three bands and are thus very unreliable (EIRS-19,40,52,55). 

Redshifts from optical spectroscopy are also available for 37 
ELAIS-IRS sourc es, most of them from s pectrosco pic follow-ups 
of ELAIS QSOs jAfonso-Luis et al.|[2004 : La Fran ca et al.ll2004l) 
or the Sloan Digital Sky Survey (SDSS? York et alia OOO*). 

Since an optical photometric or spectroscopic redshift is re- 
quired for source selection, the sample is biased against the most 
obscured sources; and the lack of constraints in the IR-to-optical 
flux ratio provides a large number of unobscured AGN which are 
absent in other high-z MIR-selected samples (see 96. lb . 

39 sources are classified as QSOs by their optical spectral en- 
ergy distribution (SED) or spectrum, while 28 are galaxies (includ- 
ing starbursts and obscured AGN) and the remaining 3 are of un- 
determined type. QSOs and galaxies can be easily distinguished by 
their optical-to-MIR flux ratios (Fig.[T]( with little overlap between 
populations. 

In order to provide a good reference sample for redshift esti- 
mation and calibrate the MIR diagnostic tools, we have compiled 
a set of MIR spectra and MIR-to-FIR photometry containing 137 
well-known low-redshift sources (the 'Library'), including QSOs, 
Seyfert I and Seyfert 2 galaxies, ULIRGs and starbursts. The MIR 
spectra are high- and low-resolution IRS spectra published in the 
literature (see Tabled, while redshifts, classification and IR pho- 
tometry for the sources were extracted from the NASA Extragalac- 
tic Database (NED). 



3 OBSERVATIONS AND DATA REDUCTION 

All targets were observed in the two low-resolution modules of the 
IRS, Short-low (SL) and Long-low (LL), exposing the two sub-slits 
on each of them (SLl, SL2, LLI and LL2) and thus covering the 5- 
38 ^im range. Two nod positions were observed in each configura- 
tion in order to provide optimal background subtraction. Exposure 
times for the LL were defined using ISOCAM 15-|,im flux densities 
in order to assure an adequate S/N at this wavelength. Integration 
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Figure 2. Residuals obtained after subtracting the broad-band IR photom- 
etry to the synthetic photometry from the IRS spectrum for each of the 
calibration bands; 8 |J,m (top), 15 |J,m (centre) and 24 |j,m (bottom). 



times for the SL were calculated from the 15-|J,m flux density as- 
suming an M82-like SED. The observations log is shown in Table 
□ 

Raw IRS data were reduced by the IRS pipeline (ver. S 15.3.0) 
at the Spitzer Science Center. This included ramp fitting, dark sub- 
traction, flat-fielding and wavelength and flux calibration. 

Since our sources are faint, additional reductions were car- 
ried out starting from the Basic Calibrated Data (BCD) stage. 
Images with the object centred in alternate sub-slits of the same 
module were subtracted (i.e: SL1-SL2, LLI-LL2) to remove sky 
background, and their corresponding uncertainty and mask im- 
ages were, respectively, added and bit-wise XOR combined!] One- 
dimensional spectra for both nod positions were extracted using the 
Spitzer IRS Custom Extraction (SPICE) version 1.4.1. The extrac- 
tion aperture was set to the SPICE default, and the Optimal Cali- 
bration Mode, in which the signal in each pixel is weighted with its 
uncertainty, was selected. 

Individual spectra were flux-calibrated with SPICE using stan- 
dard calibration files for pipeline S 15.3.0, and then combined into 
a single spectrum using custom IDL routines. 

Comparison with IRAC 8.0-[,im and MIPS 24-|,im photometry 
from SWIRE and ISOCAM 15-nm from ELAIS indicates that the 
relative flux calibration uncertainty is around 14, 35 and 8 per cent 
at 8, 15 and 24 \im, respectively. Note that the uncertainty at 15 
\im increases sharply around 1 mJy as it approximates the detec- 
tion limit of the ELAIS survey (Fig.|2j. The 8-^im fluxes from IRS 
data are in average 10 per cent lower than in the IRAC photometry, 
but we have preferred not to make any flux correction to the IRS 
spectra since the 8-[xm band covers a somewhat noisy region of the 
Short-Low spectrum in the overlapping wings of the SLl and SL2 
transmission profiles. 

Wavelength calibration is performed by the IRS pipeline, and 
is accurate to ~l/5 of the resolution element or 6A/A ~ 0.001. Re- 
duced and calibrated spectra for each of the 70 ELAIS-IRS sources 
are shown in Fig.|3] 



In the mask images, individual bits are set for each pixel that meets a 
particular condition. In the combined mask, a bit is set if it is set in at least 
one of the two original masks. 
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Figure 3. Observed IRS spectra for the ELAIS-IRS sources. The shaded ai'ea represents the lo" uncertainty in fy, while the solid line represents the smoothed 
spectrum after a Gaussian filter with a FWHM of 3 resolution elements has been applied. The error bars mark the flux density measured in the IRAC (3.6, 
4.5, 5.8 and 8 |.im) and MIPS (24 |.im) bands from SWIRE and ISOCAM (15 ^m) from ELAIS. The IR classification obtained in g]2]is indicated as QSO 
(unobscured AGN-dominated), AGN2 (obscured AGN-dominated) or SB (starburst-dominated). A question mark besides the redshift estimate indicates the 
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Figure 3. Continued. 



6 A. Herndn-Caballero et al. 




A observed (/im) 



Figure 3. Continued. 
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Figure 3. Continued. 
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Figure 3. Continued. 



4 RESULTS 

4.1 Redshifts from the IRS spectra 

Spectroscopic redshifts (z^pcc) are available for roughly half of the 
sample, all of them optical QSOs except for one obscured AGN 
(EIRS-7). The rest of the sources were selected based on photo- 
met ric redshifts (Zphot) published in the ELAIS bandmerged cata- 
log jRowan-Robinson et alj|2004l) . which were later revised with 
the addition of IR photometry fro m SWIRE and UKIDSS after 
the E LAIS-IRS sample was selected jRowan-Robinson et alj2004 
I200S) . leading to revised z < 1 in some cases. 

In many of the sources in the sample a redshift estimate can 
also be obtained from the IRS spectrum based on the PAH features, 
the silicate profile and the shape of the continuum SED (zirs)- To 
this aim, we have developed an algorithm that compares every spec- 
trum with those in the Library and determines the best-z estimate 
by a -minimization procedure (Heman-Caballero et al. in prep.). 

Reliable zirs are obtained for 26 sources, while for another 
28 a z estimate is obtained, but it is uncertain due to poor SED fit- 
ting or multiple solutions. Most sources with reliable zirs are galax- 
ies, because their PAH-emission and/or silicate-absorption features 
produce a single sharp spike in x^i^)' while in QSOs a smoother 
spectrum allows for multiple, weaker solutions. 

A comparison of optical- and MIR-spectroscopic versus pho- 
tometric redshifts (Fig. O indicates that Zphot estimates have sig- 
nificant uncertainties, larger than those of the general ELAIS sam- 

i ile but similar to the uncerta inty found in the z > 1 subsample 
Rowan-Robinson et al. Photometric redshifts are signifi- 

cantly less accurate in optical QSOs than in galaxies, due to aliasing 
problems, and may be also affected by variability issues, because 
the optical photom etry from WES was not obtained in the same 
epoch for all bands lAfonso-Luis et alj2004h . 

Since zirs estimates are stronger in galaxies and obscured 
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Figure 4. Compaiison of the photometric redshifts with those obtained from 
optical or MIR spectroscopy for sources with reliable redshift estimations. 
Triangles represent sources classified as QSOs in the optical, while squares 
represent galaxies. The open square represents EIRS-7, the only optical 
galaxy with redshift from optical spectroscopy. The dotted lines enclose 
photometric redshifts with errors lower than 10 per cent in 1 + z. 



AGN, while Zspec is more appropriate for QSOs, the combination 
of both techniques allows us to obtain reliable redshifts for 57 of 
70 sources (Table[5). The remaining 13 rely on Zphoi (EIRS-47, 60 
and 65) or unreliable zirs estimates (10 sources), and will be dis- 
carded in all redshift-sensitive analysis. 

The redshift distribution of the ELAIS-IRS sample is shown in 
Fig.H] Starburst galaxies concentrate around z ~ 1 because the se- 
lection at 15 [xm favours sources with a strong 7.7-[xm PAH feature 
at this redshift, while simultaneously prevents selection of higher- 
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Figure 5. Redshift distribution of ELAIS-IRS sources as a function of IR 
class: (a) unobscured AGN, (b) obscured AGN, (c) starburst galaxies, (d) 
all sources. The stripped histograms represent both reliable and unreliable 
redshifts, while the solid ones include only unreliable redshifts from Zphot 
orziRs. 



redshift sources since a typical starburst SED dwindles quickly for 
A <6 [im. This behaviour is also observed in 24-|xm sel ected star- 
burst galaxies, which tend to concentrate around z ~ 1.7 jYan et al] 
l2007l : lFarrah et alj 12008*). The redshift distribution for QSOs is 
smoother, with most of them in the 1 < z < 2 range, but the number 
of obscured AGN decreases sharply at z > 1.2 because of their red 
SED shortwards of 6 \im restframe. 



4.2 Monochromatic fluxes and luminosities 

We calculate restframe monochromatic fluxes (fy) and luminosities 
vLy(/i) at several wavelengths (2.2, 5.5, 7, 10 and 15 |.im). The pro- 
cedure used depends on the observed and restframe wavelength: for 
the 2.2 |.im (restframe K band) we interpolate the photometry in the 
IRAC bands (useful for sources at 0.6 < z < 2.65); for 5.5, 7 and 10 
|xm the spectrum is averaged over a narrow range (5.3-5.8, 6.6-7.4 
and 9.5-10.5 (xm respectively); while for the 15-[xm luminosity the 
procedure is more complex: in z < 1.4 sources it is obtained by av- 
eraging the 14—16 i-im spectrum; in 1.4 < z < 2 sources 15 ^m rest- 
frame is outside the spectral range observed by IRS, but not too far, 
so the 15-[xm luminosity is obtained by extrapolation of the model 
spectrum fitted in the spectral decomposition analysis (see g5.lt . In 
z > 2 sources the 15-|.im luminosity is interpolated between the red 
end of the IRS spectrum and the MIPS70 photometric point (see 
Fig.|6]for an example of the 15-|,im flux measurement in the last 2 
cases). Table|6]summarizes the results obtained. 



4.3 PAHs and Silicate features 

In the literature we find two main methods for measuring the flux 
in the aromatic bands, which differ in the way the continuum and 
the PAH features are measured: a) the interpolation method, in 
which the continuum und erlying the aromatic band s is estimated 
by linear (e.g.jR igopoulou et al.lll999h or spline (e.g. lVermeii et al.l 
I2OO2I ; ISpoon et alj,2007.) interpolation between PAH-free regions 
of the spectrum, and the flux over the estimated continuum is at- 
tributed to the PAH feature; and b) the Lorentzian method, in which 
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Figure 6. Examples of the procedure followed to estimate the restframe 
15-|im flux for sources at 1.4 < z < 2, in which the model fitted in the 
spectral decomposition is extrapolated up to 16 |-im restframe (top); and 
z > 2 sources, for which a linear interpolation between the end of the IRS 
spectrum and the MIPS70 photometric point is used {bottom). 



a Lorentzia n or Pru de profile i s assumed for the PAH features 
jLaurent et al. 2000; Smith et al . 12007; Sajina et al. 2007). 

We implemented the later method using a procedure in which 
a double power-law for the continuum and Lorentzian profiles for 
the PAH features are fitted si multaneously. The fitti ng algorithm 
is similar to that described in lSmith et al.l ( l2007h and lSaiina et al.l 
f2007) and is implemented as follows: the continuum emission is 
modeled by a linear combination of two power-laws with spectral 
indices a, p ranging from -5 to 5. The PAH features are modeled 
by Lorentzian profiles L(b,a),Ao) for the five main features (6.2, 7.7, 
8.6, 11.3 and 12.7 \im), while the rest are ignored. The amplitude 
of the Lorentzian (b) must be non-negative, its width (ai) is con- 
strained between 0.5 and 2 times the typical value found in the 
literatur^ and the central wavelength (Aq) can fluctuate up to 5 per 
cent from its expected value. 

Extinction in the continuum is implemented assuming a screen 
model, with opacity t(A) following the Galactic Centre extinc- 
tion law (GC; IChiar & Tielenslliooel) . which has been widely used 
in ot her ULIRG samples w ith moderate or high sil icate absorp- 
tion jporster Schreiber et al . 2003; Sajina et al. 2007; Polletta et al] 
Thus we fit the restframe spectra in the 5-15 [im range to 
the analytic expression: 

IV 

FJA) = {AA" + BAl')e-"^'^ + ^ L(fo„ w„ Ad{A) (I) 
We quantify the intensity of the silicate feature using the 



^01 = 0.2, 0.7, 0.4, 0.2 and 0.3 [im for the feature s centred at 6.2, 7.7, 8.6, 
11.3 and 12.7 |.im, respectively iu & Drainel200ll) . 
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'silicate strength', S ^ji (e.g. IShi et alj |2006| ; ISpoon et al] l2007l : 
iMaiolino et allboOTl) , defined as: 



(2) 



where C(Aq) = AA" + BA^ is the extinction-corrected contin- 
uum at 9.7 |.im, and F,i(/io) is the observed flux (averaged in a nar- 
row interval to minimize the effects of the noise in the spectrum). 
Ssii is positive for sources with silicates in emission and negative in 
absorption. 

In two high redshift sources (EIRS-30 and EIRS-63) 5sii could 
not be measured because 9.7 |.im restframe is outside of the ob- 
served wavelength range, and in three starburst sources (EIRS-13, 
25 and 32) the mean flux in the 9.5-9.9 |im range is negative in 
their noisy spectra. For these five sources we assume Ssii = -Tgj, 
where Tg 7 = fT(/i=9.7|xm) is the estimated apparent optical depth at 
9.7 ^im from Eq.[T] 

In most sources with silicates in emission, the silicate fea- 
ture is wider and centred at longer wavelengths (~10.5 |,im) 
than it is found when the feature appear s in absorption (see 
ISiebermiorgen et aill2005l : lNetzer et al.ll2007l for tentative explana- 
tions), and the extra flux in the silicate feature tends to increase the 
continuum level in the fit. These effects lead to a systematical un- 
der estimation of the silicate strength in sources with silicates in 
emission. To compensate for this, we run a second fitting loop in 
sources with T9 7 ~ and 5sii > in which the 9.5-1 1 |.im range is 
excluded from the fit to avoid overestimation of the continuum, and 
Ssii is measured at 10.5 |.im. 

Table |7] summarizes the results obtained from the fitting. 
Fluxes of PAH features with S/N below Icr are considered as non- 
detections and shown as upper limits. 




10 100 
X restframe [/im] 



1000 



Figure 7. SEDs models used to estimate IR luminosity of the ELAIS-IRS 
sour ces. The templates S trong FIR QSO, Weak FIR QSO and no-FIR QSO 
from lNetzer et alj )2007l) have been extended longwards of 70 |.im using the 
QSO-high template from lPolletta et al.l i2007h . 



4.4 Infrared Luminosity 

IR luminosity (Ljr) is estimated by integrating in the 8~I000 |,im 
range the template SED that best fits the MIR and far-infrared 
(FIR) photometry of each source. Template SEDs are selected from 
well-known theoretical and semi-empirical models in the litera- 
ture, including four FIR models (torus, cirrus, M82 and Arp220) 
fromH owan- Robinsod {20(31); 

one early-type, seven late-types, 
three starbursts, six AGN and three SB-l-AGN composites from 
IPoUetta et al] ( l2007h : a s et of 64 theore tical galaxy SEDs (from 
quiescent to active) from iDale & Heloul (2002); theoretical mod- 
els of normal and s tarburst galaxies as a function of Ljr from 
iLagache et alj j2003i i2004); two empirical mean SEDs (radio -loud 
and radio-quiet) from a sample of 47 QSOs telvis etaljl994h ;and 
three mean SEDs of optically-selected QSOs with strong, weak or 
no detection in the FIR (Netzer et al. 2007). 

To reduce the number of template SEDs, the best-fitting tem- 
plate for each source in the Library was obtained by fitting its 
mid-to-far IR photometry (12-850 ^irn). Template SEDs that pro- 
vided the best-fitting solution for less than two sources were dis- 
carded, and fr om very similar ones (e.g. Dale25 and Dale26 from 
bale & Heloij[2(j02l) only one of them was selected. This reduced 
the s ample of templates t o 14 m odels (see Fi g. [Til: M82 and Arp220 
from iRowan- Robinsod J200lh: GLI2 from iLagache et all j2003l. 
l2004h ; Dale26 from iDale & Helod {2OO?); NGC6240, Mrk231, 
119254s, Sey fert2, QSO-norm, QSO-high and QSO-low from 
IPoUetta et alj(l2007h: Strong F IR QSO, Weak FIR QSO and no-FIR 
QSO from dNetzer et alj|2007h . 

The best-fitting template for the ELAIS-IRS sources is deter- 
mined by fitting the IRAC 3.6, 4.5, 5.8 and 8.0 ^m, ISOCAM 15 



1,1m, and MIPS 24, 70 and 160 \im photometry, as well as synthetic 
photometry in the 10-38 f^im range obtained from the IRS spectrum. 
The synthetic filters have Gaussian transmission profiles centred at 
A, = 10.4, 12.4, 14.9, 17.9, 21.5, 25.8 and 30.9 (im, with FWHM 
= AJ6. The photometric bands shortwards of 3.6 \im are discarded 
because at z > 1 the flux in these bands is dominated by the stellar 
population in starbursts and obscured AGN, or by the accretion disc 
in QSOs, and is poorly correlated to the FIR emission from dust. 
Fitting results are shown in Fig. [8] The fits look reasonable for most 
sources, except for a significant under estimation of the MIR emis- 
sion in some obscured AGN with a steep slope in the IRAC bands 
(e.g. FIRS- 12, EIRS-47, EIRS-55). These sources fit an M82 tem- 
plate, probably because of the lack of a suitable obscured-AGN 
template. 

If we assume that the FIR SED of the ELAIS-IRS sources 
is not very different to that found in the sources in the Library, 
the uncertainties in Ljr for ELAIS-IRS can be estimated using the 
following procedure: 

1) calculate, for the sources in the Library, the best-fitting SED and 
its associated IR luminosity (Ljj^), using the IR photometry up to 
15, 25 or 60 |,im (which roughly corresponds to 38, 70 and 160 f^im 
for sources at z ~ 1.5, see Fig.|9}. 

2) compare ii^to the 'true' IR luminosity (Ljr) obtained by fitting 
all the FIR and sub-mm photometry available. 

The mean value of the logarithmic error in Lir, 
Alog(LiR) = |log(LiR) - log(L[j()|, is around 0.05 when LJ'j.is 
estimated using all the available photometry up to 60 [xm (equiv- 
alent to detection up to 160 [xm in a z ~ 1.5 ELAIS-IRS source), 
0.15 when using the photometry up to 25 |,im (equivalent to 70 [xm 
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Figure 8. Full optical to FIR SED of the ELAIS-IRS sources. The broken line represents the IRS spectrum smoothed using a boxcar filter 5 resolution elements 
wide, while eiTor bars and aiTows represent photometric measurements in the optical (I/, g, r, i, Z), NIR (J, K), IRAC (3.6, 4.5, 5.8, 8.0 \im), ISOCAM (15 
|im) and MIPS (24, 70, 160 |.im). The smooth continuous line represents the best fit template SED for the IR photometry. 
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Figure 8. Continued. 
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Figure 8. Continued. 
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Figure 9. Correspondence between photometric bands used for IR lumi- 
nosity estimations in the Library and ELAIS-IRS sources. For ELAIS-IRS 
sources aX z ~ 1-5 the MIPS 70 and MIPS 160 bands roughly correspond to 
IRAS 25 and 60 i-im photometry in local ULIRGs. The restframe spectral 
range common for both samples in IRS spectroscopy (5-15 |.im) is high- 
lighted. 



detected ELAIS-IRS sources), and 0.23 if only photometry up to 
15 |J,m is used (equivalent to ELAIS-IRS sources undetected at 70 
and 160 \.im). 

In addition to these errors, uncertainty in the MIPS70 and 
MIPS 160 photometry (up to 50 per cent in sources near the de- 
tection limit) and the bias introduced by the use of a discrete set of 
templates must also be taken into account. Considering these extra 
error sources, a factor 2-3 of uncertainty seems reasonable for most 



Figure 10. Histogram of LiRfor the ELAIS-IRS sources as calculated from 
SED fitting. The dot-filled histogram represents the whole sample, while 
those of the subsamples with MIPS70 and MIPS 160 detection are drawn in 
diagonal and horizontal lines, respectively. 



sources, significantly lower than the factor 5-10 estimated for ex - 
trapolations using only the MIPS24 photometry ( jPale et al.ll2005l) . 

Most ELAIS-IRS sources have estimated IR luminosities in 
the ULIRG range, with some extreme cases in the Hyper-luminous 
range (Lir> 10'^ Lq) and 1/4 of the sample in that of LIRGs. Fig. 
1 101 shows a histogram of the luminosity distribution of the sample. 
It is noteworthy that the sources detected at 70 and 160 |,im con- 
centrate in the higher luminosity tail of the distribution. This is an 
expected outcome given the relatively shallow depth of the MIPS70 
and MIPS 160 observations, and as a consequence, the most accu- 
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rate Ljr estimations are obtained in tlie most luminous sources. Ta- 
ble [13] indicates the best-fitting SED and estimated IR luminosity 
for the ELAIS-IRS sources. 



4.5 Star Formation Rates 

For sources dominated by star formation, the most robust mea- 
surement of the star formation rate is obtained from the integrated 
8-100 |,im IR luminosity (Lir). For starburst galaxies, iKennicuta 
( Il998h finds the relation: 



1.72 X 10-'"LiR [Lq] 



(3) 



Unfortunately, Lir has big uncertainties for individual sources, 
and in sources with AGN activity there is not a straightforward 
method to isolate the starburst and AGN contributions to Lir. 

Many recent papers find a strong correlation between the 
luminosity of the PAH bands (Lrah) and Lir in starburst- 
dominated sources (e.g . iRigopoulou et'^ 1 19991 : | Peeters et al.l 
l2004l :l Brandl et al ] |2006n . The ratio Lpah/^ir depends on the prop- 
erties of the interstellar medium, and thus it varies with the lumi- 
nosity of the sources: in normal spiral galaxies, with quiescent star 
formation, the mean ratio is highest (Smith et al...2007.) . in starburst 
galaxies it is somewhat lo wer, and in local ULIRGs with massive 
star formation, it i s lowest jRigoDoulouet^l 19991 : ILu et al.ll2003l : 
iNetzer et al.ll2007l) . Nevertheless, the dependency of LpAn/iiR with 
Lir is small r elative to the dispersion between sources of the same 
luminosity: iRigopoulou et al.l ( Il999i) find for the 7.7-|^im PAH a 
mean value of 8.1 x 10"^ in starburst galaxies a nd 5.5 x 10"^ in 
starburst-dominated ULIRGs. lSrandl et all ^2006^ find that the 6.2- 
[xm PAH luminosity relates to Lir trough a power-law of index 1.13 
in a sample spanning 2 orders of magnitude in Lir. 

Some authors use the 6.2-|xm PAH because its measurement 
suffers little contaminatio n from other MIR features jPeeters et al.l 
l2004l : iBrandl et alj [2006^ while others prefe r the 7.7-M,rn PAH 
because it is th e strongest PAH band (Ri gopoulou et al.l Il999l : 
iLutz et al. I l2003h . or even use the combined flux of two or more 
features (e.g. lLu et alj2003l : lFarrah et al.ll2008l) . 

The estimated luminosity in the PAH features dep ends 
strongly on the measurement procedure: ISmith et^ ( l2007h find 
discrepancies of a factor 1.7 at 6.2 \\.m. and 3.5 at 7.7 [im depend- 
ing on whether the continuum is calculated by spectral decompo- 
sition into Lorentzian-Fcontinuum or it is interpolated between ad- 
jacent regions. This suggests that the discrepancies in the value of 
^pah/^ir found by different authors (see Table 1 121) are mainly due 
to systematics. In addition, the integrated infrared luminosity is not 
the same throughout th e literature: some a uthors define Lir in the 
range 1-1000 \.\.m (e.g. 'Farrah et al !"2008^, while ot hers use only 
the FI R (42-122 |xm, Lu et al. (2003); 40-500 um. IPeeters et al.l 
( l2004l) ). 

Because of these issues, we obtain our own estimate of the 
average Lpah/^ir in the starburst-dominated ELAIS-IRS sources 
before using the PAH luminosity to estimate SFR in the whole 
sample. The mean values of Lrah/^ir found in the 7 starburst- 
dominated sources with no sign of significant AGN emission 
(EIRS-21, 25, 37 and 41 excluded) are: L62/L1R = 0.012 + 0.007, 
L77/LiR= 0.038 ± 0.016 and Lhs/Lir^ 0.011 ± 0.007. 

The ratio found for th e 7.7-|xm PAH is similar to that ob- 
tained by L utz_gt al.) j 2003l) in a sample of starburst galaxies and 
ISmith et all 12007) in nor mal and starbu rst galaxies. It agrees with 
that found by IRigopoulou et alj ( Il999h in starburst galaxies and 
ULIRG s if we take into a ccount the factor 3.5 decrement calcu- 
lated bv lSmith et al.l ( |2007|) for fluxes measured by interpolation of 
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Figure 11. Histogram of star formation rates as estimated from the lumi- 
nosity of the 6.2, 7.7 and 11.3 \im PAH features in individual ELAIS-IRS 
sources. The solid histogram represent SFR estimates in sources where at 
least two PAH features are detected, while the stripped histogram indicates 
upper limits for sources with detection in just one or none of the PAH bands. 



the continuum using splines. For the 6.2-|xm feature, our ratio is in 
agreement with those of Smith et al. ( 2007) and Farrah et al. ( 200^ 
but it is roughly double the ratio found by Spoo n et alj J200 4) in a 
sample of normal and starburst galaxies even after correction for 
the differences in the continuum estimation procedure. 

By applying the calibration in Eq. [3] and assuming that the 
properties of the interstellar medium in which the PAH emission 
originates are not significantly altered by the AGN activity, the SFR 
in the whole sample can be calculated as: 



5Fi?[Meyr-'] 
5Fi?[MQyr-'] 
5F^[Moyr-'] 



1.40X 10-«Lf,2 [Lo] 
4.56 X 10-" L77 [Lq] 
1.52 X 10-' Li,3 [Lq] 



(4) 
(5) 
(6) 



For sources in which two or three PAH features are detected 
with S/N > 2 we estimate SFR from their combined luminosity, 
while if only one or no PAH features are detected, we average esti- 
mates from individual upper limits. 

Table [13] shows the estimated SFR for ELAIS-IRS sources. 
They range from ~10 to ~1000 Meyr"'. The stated uncertainties 
account for the uncertainty in the PAH flux only; if we consider the 
uncertainty in the IR luminosity and the dispersion in Lpah/Lir of 
the starburst-dominated sources, a more realistic estimate would be 
a factor 2 uncertainty in SFR. 

Figure [TT] shows de distribution of SFR in the ELAIS-IRS 
sample. Most estimates are in the range 150-600 Moyr"', but if 
we consider upper limits for sources with undetected PAHs, which 
account for more than half of the AGN population, many of them 
have probably SFR < 100 Moyr"'. 



5 SOURCE CLASSIFICATION 
5.1 Spectral Decomposition 

A powerful tool for diagnostics is decomposition into several spec- 
tral components, which c an provide consi derable insight to the 
physics of the sources (e.g. lTranetal.l200ll) . 

MIR emission from active galaxies arises mostl y from HII re- 
gions, Photodissociation Regions (PDRs) and AGN jLaurent et al.l 
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Figure 12. Relative contribution to the integrated 5-15 restframe lu- 
minosity of the AGN and HII spectral components as derived from the 
spectral decomposition for sources in the Library and ELAIS-IRS samples. 
Sources from the Library are plotted with different symbols according to 
their optical-spectroscopic classification from NED: QSOs and Seyfert 1 
galaxies (triangles), obscured quasars and Seyfert 2 galaxies (inveiled tri- 
angles), starburst galaxies (circles), and ULIRGs (squares). Solid symbols 
indicate sources with a strong starburst component, while open symbols 
represent mostly pure-AGN sources. Symbols for the ELAIS-IRS sources 
indicate their MIR classification: unobscured AGN (plus signs), obscured 
AGN (diamonds) or starbursts (asterisks). Dotted lines indicate the relative 
contribution of the PDR spectral component to the 5-15 |.im restframe lu- 
minosity. 



l200d) . so we have used a simple model comprising the superposi- 
tion of three spectral templates (AGN, PDR and HII) obscured by 
a screen of dust. 

The spectral decomposition is performed by fitting the 5-15 
[xm restframe range of the spectrum to a parametrized F((/i) of the 
form: 



""""'(fli/AGNW + fla/raiW + a3/pDR(^)) 



(7) 



where the free parameters (b, ai, a2, aj) are calculated using a 
Levenberg-Marquardt ^'-minimization algorithm, and must have 
non-negative values. 

To keep the number of free parameters low, a single optical 
depth value is applied to the three spectral components, but for most 
sources the solution is almost identical when two optical depth val- 
ues (one for the AGN and other for the PDR and HII components) 
are used, with an increased value due to the extra free parameter. 

t(A) is obtained from the GC extinction law, while /agn is 
repres ented by the IRS spe ctrum of the Seyfert 1 galaxy NGC 
3515 jBuchanan et"ai]|2006l) , /pdr by the ISOCAM spectrum of 
a PDR in the reflection nebula NGC 7023 jCesarskv et alJI 19963) 
and /hi: by the ISOCAM sp ectrum of M17 in the vicinity of OB 
stars dCesarskv et aljl996bh . 

We quantify the contribution of each spectral component to 
the MIR spectrum as the ratio of its integrated luminosity to total 
luminosity in the 5-15 ^m restframe range, ^agn, ''hii and tpdr for 
the AGN, HII and PDR components. Figure [12] shows the distri- 
bution of spectral components for the Library and the ELAIS-IRS 
sample. 



5.2 Diagnostic Diagrams 

A number of MIR diagnostic diagrams have be en proposed 
to determine the nature of infrared galaxies (e.g. Genzel et al.l 



1998VLutz et al." 1998V Rigopoulou et al.ll999l : lLaurent et al.l2OO0l: 
Peelers et al. 2004; Spoon et al. 2007). They rely mainly on the 
strength of the PAH features to infer the dominant energy source 
(AGN or starburst), and on the slope of the continuum and sign 
of the silicate strength (emission or absorption) to distinguish ob- 
scured and unobscured sources. 

There is a strong correlation between the optical and MIR 
spectrum of AGN: type 2 AGN usually show a red MIR spectrum, 
and the silicate feature in absorption, while most type 1 AGN show 
a bluer MIR spectrum, with stronger NIR emission, and the silicate 
feature absent or in emission. Nevertheless, the match between pop- 
ulations is not perfect, and many mismatches (most notably, type 1 
AGN with silicates in absorptio n and type 2 with silicates in emis- 
sion) have been reported (e.g. [Buchanan et al]|2006l : iBrand et al.l 
2007; Haoet al. 2007). 

A silicate emission feature in type 2 AGN can be repr oduced 
by a low-Si, low-tau smooth torus model jpritz et al]|2006l) . or by 
silicate emission being pro duced in an extend ed region, such as a 
dusty Narrow Line Region dSturm et alj|2005h . Silicate absorption 
in type 1 AGN can be observed if the Broad Line Region is partially 
obscured, m aybe by dust in the galaxy. A very remarka ble case 
is Mrk 231 jLonsdale et ai]|2003bl; IWeedman et al.ll2005h , which 
shows a br oad emission-line spectrum reddened with Ay ~ 2 in 
the optical jBoksenbe rg et al. 1977; Lipari et al.' '1994"), and strong 
silicate absorption in the MIR (S poon et al, 200 2). 

Since making assumptions on the Seyfert type of AGN based 
on their IR properties alone would be misleading, we will refrain 
from doing that and will instead classify the AGN as obscured or 
unobscured in the IR. 

Using some well known diagnostics along with some others 
from our own we separate the MIR spectra of the ELAIS-IRS 
sources into three main categories: starburst-dominated (SB), un- 
obscured AGN-dominated (AGNl) and obscured AGN-dominated 
(AGN2). Every diagnostic criterion is calibrated using the known 
nature of sources in the low-redshift Library, including Seyfert type 
for AGN, and a final MIR classification is obtained as the most re- 
peated result among all criteria. 

Measurements on the ELAIS-IRS and Library spectra were 
obtained using exactly the same procedures, to prevent distortions 
in the diagnostic diagrams due to systematics. Limits between pop- 
ulations in the diagnostic diagrams are determined based solely on 
the Library spectra, and are fine-tuned to minimize the number of 
misclassifications within the Library. 

Figure [13] sh ows an adapta t ion o f three well known diag- 
nostic diagrams by [Laurent et al.l feOOOl) . [Peeters et al.[ ( [2004[) and 
ISpoon et al.[ ( l2007h . as well as a new one that we propose. 

In the diagnostic diagram o f [Laurent et al] ( |2000|) the 6.2-|im 
PAH equivalent width (EW(,2) easily separates starburst-dominated 
sources from the AGN-dominated ones, while the ratio between 
15 and 5.5 [.im restframe spectral luminosities (L15/L55) distin- 
guishes sources dominated by hot dust from an unobscured AGN 
(low L15/L55) or cold dust by an obscured AGN or starburst (high 
Lis/Z-ss). The value of £1^62 that best separates starburst-dominated 
and AGN-dominated sources from the Library is ~ 0.2, while in 
L15/L55 we find too much overlap between populations and prefer 
to drop it as a diagnostics criterion. 

In the diagram from [Peeters et al.l ( [2004[) the ratio between L55 
and LiR is related to the relative amount of hot and cold dust in the 
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Figure 13 . MIR diagno stic diagrams applied to sources in the Library and ELAIS-IRS, adapted from lLaurent et al] feOOOl) (top fgffl. lPeeters et al.^ 1 2004) (top 
right) and lSpoon et alj feoO?.) (bottom left). A new diagnostic diagram based on the 7.7-|im PAH line-to-continuum ratio and 15 |J,m to 10 |J,m luminosity ratio 
is also shown (bottom right). Symbols are the same as in Fig. [12] Dotted lines represent limiting values between populations for parameters used as diagnostic 
criteria (see Table[8). Dashed lines in the Laurent diagram represent mixing lines between the AGN (NGC 3515), HII (M17) and FDR (NGC 7023) templates 
used in the spectral decomposition discussed in jl5.ll 



galaxy, being higher for unobscured AGN and lower for starbursts 
and deeply embedded sources. The 6.2-|xm PAH luminosity (L52) 
to Lir ratio depends on the relative contribution of unobscured star 
formation to the IR luminosity, and is higher for unobscured star- 
bursts or AGN+SB composite sources, while lower for the dust 
enshrouded ones. Most ELAIS-IRS sources have too high upper 
limits on this ratio to be meaningful, so we also drop this parameter 
as diagnostic criterion. 

The diagram from lSpoon et is based on EW^2 and 

the strength of the silicate feature. Galaxies populate their di- 
agram along two branches: one horizontal branch from AGN- 
dominated to starburst-dominated spectra, and one diagonal branch 
from highly obscured nuclei to starburst-dominated spectra. They 
hypothesize that the separation into two branches reflects two dif- 
ferent configurations of the obscuring dust: clumpy in the horizon- 
tal and non-clumpy (homogeneous shell) in the diagonal. In our 
adaptation of the Spoon diagram (Fig.[T3] bottom left) it is notewor- 
thy the lack of ELAIS-IRS sources in the high-extinction, weak- 
PAH region of the diagram, where many local ULIRGs are located. 
The silicate strength is used to separate obscured and unobscured 
AGN, but most sources have a small absolute value of Ssh, and its 
uncertainty can lead to a wrong classification in some cases. 



We propose a new diagnostic diagram, inspired by the one 
from Spoon but optimized for low S/N sources (Fig. \T3\ bot- 
tom right). In our modified Spoon diagram, EW(,2 is replaced by 
the line (peak) to continuum ratio of the 7.7-[im PAH feature 
(l/c-jj). This parameter h as been extensively used to distinguish star- 
bursts from AGN (e g. iLutz et al.l [1993 : iRigopoulo u et al.lll999l : 
IWeedman et"ai]|2006h . and works better in low S/N sources be- 
cause the 7.7-|xm feature is stronger than the one at 6.2 [xm. The 
silicate strength is replaced by a more straight-forward and model- 
independent measurement: the ratio between Lis and the spectral 
luminosity averaged in a narrow band centred at 10 |xm restframe 
(Lio). This ratio separates type 1 and 2 AGN from the Library more 
elficiently than 5sii does, because it combines the two main distinc- 
tive factors: absorption/emission in the silicate band and slope of 
the continuum. In addition, it is less sensitive than 5sii to noise in 
the spectrum. Its main drawback is that 15 (xm restframe is unob- 
servable by IRS for sources at z > 1.5, but it can be interpolated 
using the MIPS70 photometry as we showed in 34.21 Type 1 and 
2 AGN from the Library are better separated in this diagram than 
in the original from Spoon in both axes, and only the type 1 AGN 
with a strong starburst contribution overlap significantly with the 
type 2 ones. This overlap is produced by a higher slope of the con- 
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tinuum and depleted silicate feature in type 1 AGN with a starburst 
component. 

Finally, the spectral decomposition diagram in Fig.[T2]can also 
be used for diagnostics, since starburst-dominated sources locate 
near the left border, because of strong PDR and/or HII component, 
while AGN-dominated sources concentrate in the lower-right cor- 
ner. 

From the five diagnostic diagrams described, we extract six 
criteria to classify infrared galaxies (Table[8l(. 

Criteria CI to C4 distinguish AGN from starbursts based on 
the strength of the PAH features and MIR-to-FIR luminosity ratios, 
while C5 and C6 rely on the silicate feature to separate obscured 
and unobscured AGN. 

The final class assigned to each source is the most repeated 
one among the six criteria. This combination of multiple criteria 
allows for a more robust classification than that obtained using a 
single diagnostic diagram, since it is less likely to be altered by 
statistical uncertainties. This is most significant in the lower S/N 
spectra, and also in obscured AGN, which some criteria often clas- 
sify as unobscured AGN or starbursts because typical values for 
obscured AGN lay in a narrow range between those of starbursts 
and unobscured AGN. 

For sources in the Library, the IR classification is consistent 
with the optical class for more than 90 per cent of them (see Ta- 
bleO. An analysis of the results from individual criteria shows that 
mismatches are caused by two main factors: a) strong star forma- 
tion in type 1 AGN, which leads to classification as obscured AGN 
via high values of //C77 and Lu/Liq. b) silicate emission in type 2 
AGN, which confuses the two criteria based on the silicate feature. 

In the ELAIS-IRS sample we expect a somewhat higher mis- 
match rate because the lower S/N of the spectra increases the like- 
lyhood of incorrect classification by some of the individual criteria. 
In all diagnostic diagrams the ELAIS-IRS sources are plotted with 
symbols representing their final IR class, so that mismatches be- 
tween individual criteria and the final classification can be easily 
recognized. 

Based on the criteria presented above, 41 sources are clas- 
sified as unobscured AGN (AGNl), and 17 are obscured AGN 
(AGN2). The remaining 1 1 sources are classified as starbursts (SB), 
including a few sources with strong continuum (EIRS-21, EIRS-25, 
EIRS-37), which could host an obscured AGN, and an AGNl-l-SB 
composite (EIRS-41). Table 1 101 shows the classification from indi- 
vidual criteria as well as the final IR class for the whole ELAIS-IRS 
sample. 

A comparison of the optical versus IR classification for 
ELAIS-IRS sources (Table [TTJ shows 12 mismatches: seven are 
optical galaxies with unreliable redshifts classified in the IR as un- 
obscured AGN because of weak or absent PAH features or presence 
of silicate emission. Some of them (EIRS-8, EIRS-38) present a 
rather flat slope in the IRAC bands, indicating low extinction, while 
others (EIRS-47, EIRS-50, EIRS-56, EIRS-60) show a steep slope 
and FIR emission, suggesting a high level of extinction. The other 
galaxy, EIRS-27, shows a very strong silicate emission feature, 
similar to that found in the Seyfert 2 galaxy IRAS F01475-0740 
( iBuchanan et alj2006l : lHa"o et al.l2007h . The remaining five sources 
are optical QSOs with reliable spectroscopic redshifts, which are 
classified as obscured AGN in the IR due to intense PAH emission 
and a small but negative value of S^w- In these sources, the dust in 
star-forming regions may have depleted the silicate emission fea- 
ture in the AGN spectrum enough to make it appear in absorption. 
In three of these sources (EIRS-5, EIRS-39 and EIRS-54) the cri- 
terion in L15/L10 predicts AGNI while 5sii < 0. A fourth source. 



EIRS-63, is at z = 3.094 and the silicate feature moves out of the 
range observed by IRS, so the IR diagnostic is inconclusive. Fi- 
nally, EIRS-58 has a deep silicate feature and is classified as AGN2 
by both L15/L10 and 5sii, but its optical SED is blue and shows a 
power-law spectrum in the optical and UV up to the Lyman break. 

The best fitting IR SED and MIR classification are consis- 
tent with each other for most sources. All sources classified as 
unobscured AGN fit best to one of the seven QSO SEDs, except 
EIRS-47, which depends on a photometric redshift and is obscured 
in the optical. Obscured AGN also select mainly QSO templates 
(most significantly Mrk 231, a reddened type- 1 QSO with strong 
silicate absorption), but 1/3 of them fit best M82 or 119254s (an ob- 
scured AGN + starburst composite). Starburst sources favour pure- 
starburst (GLI2 and Dale26) or composite (Seyfert 2, 119254s and 
NGC 6240) templates, except for two sources with strong MIR con- 
tinuum (EIRS-41 and EIRS-37) which fit best QSO templates. 

5.3 Averaged spectra 

We have used the spectral classification in 95.2l to calculate average 
spectra of each population, which have better S/N than the spec- 
tra of individual sources. In the average spectrum we only include 
sources with reliable redshift estimations, to prevent as much as 
possible the dilution of fine details by the dispersion in z- 

To calculate the averaged spectrum, individual spectra are 
shifted to the rest-frame, normalized at 7 |.im and interpolated in 
a common array of wavelengths uniformly distributed in \og(X). 
For each value of A we discard the highest and lowest value of 
and the remaining ones are combined using a weighted average. 

Figure [14] shows the averaged spectra for starbursts, optical 
QSOs and obscured AGN (optical galaxies classified as AGN in 
the IR) in the sample. 

The quasar spectrum can be approximated by a power-law of 
index o- = -I in the range 4—16 |xm (Fig. [TS] ), which is slightly 
steeper than the index found in th e averaged spectrum of a sample 
of local QSOs jNetzer et al.ll2007h . In addition, the silicate feature 
appears in emission in the averaged spectrum, but with a very low 
profile. These findings suggest the ELAIS-IRS quasars could be 
more obscured than average, nearby optically selected QSOs. See 
36.4| for a discussion on the connection between the silicate feature 
and the MIR spectral index. 

The slope of the spectrum increases shortwards of 3 [im to- 
wards the 1 i^im minimum of the SED, as has also been observed in 
other optically selected quasar samples (e.g. lHatziminaoglou et al.l 
,2005.) . This is interpreted within the torus model as thermal emis- 
sion from graphite grains at temperatures up to ~1500 K. 

In the averaged spectrum there are tentative detections of some 
ionic lines of neon, [Ne 11] 12.81 |,im and [Ne m] 15.56 [xm, but 
since they have low/medium excitation potentials, their emission 
can also be excited by star formation. The PAH features are weak 
but detected, with EWa ~ 0.015 and EWn ~ 0.05. In |63]we es- 
timate mean PAH luminosities and star formation rates for these 
sources. 

The averaged spectrum of the obscured AGN shows a steep 
slope at /I < 8 |,im {a = -2) but it flattens longwards of the silicate 
trough at ~I0 |,im (Fig.ll6t. The depth of the silicate feature is rel- 



and QS02s (e.g. 


Sturm et al. 2005: Weedman et al. 2006: Deo et al.l 


l2007l;lHaoetal.l 


2007). However, the requirement of optical spec- 



troscopic or photometric redshift for selection may have reduced 
the average silicate absorption by biasing the sample against the 
most obscured sources. 
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Figure 14. Composite spectrum of ELAIS-IRS sources separated based on 
their IR classification: unobscured AGN {top), obscured AGN (centre) and 
starbursts (bottom). The continuous line represents the average spectrum, 
while the shaded area represents the Icr dispersion of individual spectra. 
The dashed line in the bottom panel represents the average starburst spec- 
trum after removal of sources with a strong continuum component (EIRS- 
21,25, 37 and 41). 
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Figure 15. Composite spectrum of the ELAIS-IRS QS Os compared to the 
composite spectrum of a local sample of PG QSOs from lNetzer et alj fcOOTh 
(Netzer+07) and a sample of 12 type I mm-selected QSOs at z ~ 2 from 
iLutz et alj <2008ll (Lutz+08). A power law fy x v " with index a = 1 is 
shown for reference (dashed line). 
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Figure 16. Composite spectrum for the ELAIS-IRS AGN2 sources com- 
pared to the composite spectrum of X-ray selected, higher-luminosity 
QS02s from iPolIetta et alj f2003) in which the NIR-to-MIR spectrum is 
best fitted by a torus model (PoIIetta T) or torus + cold absorber (Polletta 
T-hC). a power-law /y oc y"" with index a = 2 is shown for reference 
(dashed lines). 

Emission lines for ionic species are very weak or absent, with 
the significant exception of [Ne ii] 12.81 |xm. The 7.7-|^im PAH fea- 
ture seems rather strong, but it could be magnified by a peak in 
the underlying continuum emission at 8 ^m caused by the onset of 
the silicate absorption feature. The other PAH features are weak, 
with only the one at 11.3 |,im clearly detected. The 6.2 \\m feature 
is missing even if it is detected in many individual spectra. This 
may be caused by the cancellation elTect produced by strong 6.0 
|xm absorption ban ds by water ice or hydrocarbons in some sources 
dSaiina et alj2007h . 

In the averaged spectrum of the starburst-dominated galaxies 
we find, in addition to the usual PAH bands at 6.2, 7.7, 8.6, 11.3 
and 12.7 |,im, a few additional, fainter bands, centred at 5.2, 5.8, 
14.2 and 17.0 [xm, which are also common in higher S/N spectra 
of low redshift galaxies. The 3.3-^im band, relatively strong in ISO 
SWS spectra of many local ULIRGs, seems absent; but the IRS 
spectra of the ELAIS-IRS starbursts are very noisy in this range. 

Some ionic lines frequently observed in local starbursts also 
appear in the average spectrum, most significantly [Ar ii] 6.99 |.im, 
[Ar in] 8.99 |.im, [Ne ii] 12.81 |,im (which overlaps with the 12.7-|xm 
PAH) and [Ne iii] 15.56 |.im. It is noteworthy the absence of rota- 
tional transitions of the H2 molecule, with the possible exception of 
S(3) 9.67 |.im. The weak continuum shortwards of 6 |xm indicates 
that the relative contribution of AGN or obscured star formation 
is small. If the four starburst-dominated sources with suspected 
or known AGN contribution are discarded, the average spectrum 
of the remaining starburst sources (EIRS-2, 3, 9, 13, 14, 16 and 
32) has a lower continuum emission at /I < 6 ^im, but remains un- 
changed at longer wavelengths while the PAH features get signifi- 
cantly boosted. 



6 DISCUSSION 

6.1 Comparison with other high-z samples 

Before Spitzer, the spectroscopic study of ULIRGs was restricted 
to a few local sources due to the low sensitivity of previous in- 
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Figure 17. Colour-colour diagram for the most repr esentative high-z 
ULIRG samples observed w ith the IRS: Weedman et alj 12006)) (x signs), 
lYanetal.lf2007l) (diamonds). lFarrah et alj)2008h (plussigns), |Polletta et alj 
i2008t) Oriangles) and ELAIS-IRS (asterisks). The dashed vertical line indi- 
cates the lower limit in 524/Sr for the sample of iHouck et alj j2005l) . The 
continuous and dotted lines represent, respectively, tracks for M82 and 
Arp220 redshifted from ^ = to z = 3. 



struments; but in the last four years, thanks to the unprecedented 
sensitivity of IRS, it has been possible to obtain MIR spectra for 
several hundred ULIRGs up to z ~ 3. 

Most of these sources were selected from wide-area MIR 
surveys, and are usually heavily obscured, because only sources 
with high MIR-to- optical flux ratios were selected fo r observa- 
tion with IRS (e.g. iHouck et alJllOOSl: lYan et al.ll2007h . In some 
cases, additional constraints were imposed to favour the selec- 
tion of starburst or AG N-domin ated sources (Weedman et al. 2006; 
iFarrah et alj2008l : iPolletta et alllOO^j) . 

By contrast, the selection criteria for the ELAIS-IRS sample 
put no constraints on the MIR-to-optical flux ratio of the sources, 
but required optical spectroscopic or photometric redshifts. This 
introduced a bias that is not easy to quantify, but certainly favored 
sources with little or no obscuration at moderate redshifts. In addi- 
tion, selection at 15 |,im instead of 24 |,im also favours lower-redshift 
sources, because the PAH features and the steep slope shortwards 
of 6 i^im in obscured AGN enter the selection band at lower z- In 
this sense, the ELAIS-IRS sample can be seen as complementary 
to other higher-z spectroscopic samples. 

The selection criteria leave their fingerprint in a colour-colour 
diagram like the one in Fig.[T7] Indeed, the ratio between the flux 
density at 24 \im and an optical band, used as selection criterion for 
all samples except ELAIS-IRS, determines to a large extent the na- 
ture of the sources selected. The most extreme cases (SnlS, > 2000) 
dominate the samples from Houck and Polletta. These sources are 
strongly obscured in the optical and show steep IR spectra, with 
strong silicate absorption, and tend to be at high redshift (z > 1.7). 

Ratios of 300 < S24/5, < 2000 correspond to the obscured 
AGN in ELAIS-IRS, the starburst galaxies from Farrah and the 
composite sources from Yan. The distribution of z goes down to 
~0.5 with a gap around z ~ 1.4 because of the silicate feature en- 
tering the MIPS24 band. 

At S24/5,- ~ 100 we find the ELAIS-IRS starburst at 0.6 < z < 
1, which are bluer in 524/5^ and 524/58 than starburst sources from 
Farrah or Yan; and at 524/5^ < 100 we find most of the ELAIS-IRS 



Figure 18. MIR luminosity as a function of redshift for some high-z ULIRG 
samples. Symbols are as in Fig. 1171 The eight ELAIS-IRS sources with 
lower 5.8-^im luminosity that depart from the main trend are starbursts, 
while most of the other sources are AGN or AGN-dominated composite 
sources. 



quasars, which are slightly or not obscured, spanning a wide range 
in redshift up to z ~ 3. 

Note that excluding the ELAIS-IRS sources, the rest are sig- 
nificantly redder in 524/5^ than would be M82 or Arp220 at the same 
redshift, indicating that these sources are more obscured than many 
local ULIRGs, and probably represent only the most extreme cases 
of the ULIRG population at high redshift. 

Since most of the sources in these samples have 524 ~ 1 mjy, 
their distribution in the L-z plane covers a very narrow strip (Fig. 
list , indicating that, except for a few sources w ith ultra-deep IRS 
spectra down to ~0.1 mJy jTepUtz et al.1l2007h . only the most lu- 
minous sources at each redshift have been probed so far. A deeper 
insight into the dominant population of LIRGs at z > 1 will be pos- 
sible once the next generation of MIR spectrographs onboard JWST 
and SPICA become available. 



6.2 Starburst ULIRGs at z ~ 1 

We have classified 10 of the 27 optically-faint ELAIS-IRS sources 
as starburst-dominated galaxies, and there is also an optical quasar 
(EIRS-41) dominated by starburst features in the MIR. 

The selection criteria required Zphot > 1 for galaxies, but zirs 
estimates are in the 0.6-1.2 range for all starburst galaxies, with (z) 
= 0.8. This is caused by uncertainties in Zphot and by the selection 
at 15 |xm, which favours sources with strong PAHs at z ~ 1, while 
at higher redshift requires a strong continuum at A< 6^im, favoring 
AGN. 

In spite of this bias towards PAH-strong sources at z ~ 1 , only 
1/3 of the 0.6 < z < 1.2 ELAIS-IRS sources are starbursts, while 
2/3 are AGN. Most significantly, if we exclude the QSOs only half 
of the galaxies in this range are starbursts, the rest being obscured 
AGN. This indicates that the population of 15-[xm sources at z ~ 1 
is comprised mostly of AGN for 5i5 > 1 mJy. The starbursts could 
still dominate the highly obscured galaxy population, since they are 
excluded from the ELAIS-IRS sample because they are too faint in 
the optical; but this is probably not the case, because the AGN also 
do minate the 24-|xm coun ts of sources brighter than 1 mJy at z ~ 
1.7 dWeedman et alj2005l) . 
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Figure 19. Composite spectra for sources with strong PAHs in four high- 
z samples: E LAIS -IRS starbursts (solid line), z ~ 1.8 sta rbursts from 
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Figure 20. Location in the Spoon diagram of the composite spectra of sev- 
eral high-z starburst-dominated samples (filled stars). The open star rep- 
resents the averaged spectrum of starburst-dominated ELAIS-IRS sources 
after removing the four sources with strong MIR continuum. Positions for 
all sources in the Library are also shown as a reference. 



Figure 21. Luminosity ratio for the 6.2, 7.7, and 11.3 |-im PAH features 
in the composite spectra of starburst-dominated sources in high-z sam- 
ples (big solid stars). The open star represents the averaged spectrum of 
starburst-dominated ELAIS-IRS sources after removing the four sources 
with strong MIR continuum, while individual starburst-dominated sources 
from ELAIS-IRS are shown as asterisks. Sources from the Library with a 
starburst component are also shown as a reference: starburst galaxies (cir- 
cles), starburst ULIRGs (squares), type 1 AGN with starburst component 
(triangles) and type 2 AGN with starburst (inverted triangles). 



In the diagram of PAH luminosity ratios (Fig. I21l l the com- 
posite spectrum of ELAIS-IRS starburst-dominated sources also 
situates among local starburst galaxies, with Lj2/iii3 substantially 
higher than in typical starburst-dominated ULIRGs. This ratio is 



relate d to the ionization state of the PAH rnolecul e s (.Draine & Li 



2001 



iRapacioU etal] l2005h ISrandl et alj l2006l : iGalliano et a" 



20081) . a high ratio indicating abundance of cationic species, while 
a low ratio indicates most PAH carriers are neutral. 

^62/^77 shows anticorrelation with the extinction 
llRigopo ulou et al" I ll999l) . and is usually lower in local U LIRGs 
than in less luminous starburst galaxies jLutz et al] 1 19981) . In the 
ELAIS-IRS average, it is somewhat higher than in most local 
ULIRGs, suggesting a lower extinction level. 

Lower ionization an extinction values than in local ULIRGs 
have also been fo und in higher-z starburst-dominated galaxies 
jFarrah et al. I l200i) . and could indicate that the star formation is 
extended over larger regions relative to their local counterparts, or, 
alternatively, that their gas-to-dust ratio is higher. 



The averaged spectrum of the ELAIS-IRS starbursts is simi- 
lar to that found in higher redshift samples selected at 24 |,im (Fig. 
I19t . Differences between them in the depth of the silicate feature 
and PAH equivalent width can be interpreted as a combination of 
a 'pure starburst' (HII-)-PDR) spectrum and an AGN spectrum in 
variable proportions. In this case, starburst galaxies from Farrah 
and ELAIS-IRS would be the less contaminated by AGN dust con- 
tinuum, while in those of Yan the continuum would dominate de 
MIR spectrum. 

The Spoon diagram for averaged spectra (Fig. |20|) indicates 
that the composite spectrum of the ELAIS-IRS starbursts is in the 
range of local starburst ULIRGs. If the starburst-dominated sources 
with strong continuum are removed from the composite, the result- 
ing E'W(,2 increases threefold and resembles that of lower luminos- 
ity local starbursts. 



6.3 Star formation in the AGN 

It has been stated that the SFR o f Seyfert galaxies corre l ates with 
the lu i mnosity of the AGN (e.g. ISchweitzer et al.ll2006l : IShi et al.1 
I2OO7I ; iMaiolino et all l2007h and that the S FR is systematically 
higher in ty pe 2 Seyferts relative to type 1 jMaiolino et al.|[l99a ; 
iBuchanan et al. 2006). 

For higher luminosity AGN there are contradictory results: 
most local QSOs are hosted by ellipti cal galaxies, with very lit- 
tle star formation ( lDunlopetal]|2003h . and only 3 per cent of 
them show signs of interaction with other galaxies jGuvon et al.l 
,2006,,) . Measurements in the [O 11] 3727 -A line indicate that the 
SFR i s very low (a fe w Moyr ') in optically selec t ed QS Os jHol 
I2OO5I : iKim et al.ll200^ . In the IR, ISchweitzer elal] ( l2006h find in 
IRS spectra of Palomar-Green QSOs that more than 30 per cent of 
the IR luminosity is due to star formation, but a later work from 
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IShi et alj lioOTh reduces the starburst contribution to ~25 per cent 
of the flux at 70 and 160 \an -which translates to ~10 per cent of 
LiR if we assume templates like Dale26 for the starburst and QSO- 
norm for the AGN component- for a larger sample of optically 
selected QSOs. Using MIR and FIR photometry from SWIRE, 
iHatziminaoglou et alj | |200§ ) estimate relative AGN and SB contri- 
butions to LiR in a sample of 70 SDSS QSOs with 70-|,im detection, 
and find the AGN contribution varies from ~ 100 to ~20 per cent of 

^-IR- 

Optically selected QSOs, ho wever, represent less th a n half 
of the total Q SO population jMartmez-Sansigre et ^ I2OO5I : 
IStem et ^120051) , the remai ning being obscured quasars (QS0 2s) 
that show up in the inf rared jCutri et al. 200lVLacv et al]|2004[) ra- 
dio dWhite et alj2003h or X-rays jNorman et aU2003) . 

Conversely to what is found in local QSOs, HST images of 
a sample of 0.3 < z < 0.8 QS02s show sig ns of perturbation s 
and intense star formation in the host galaxy jLacv et al]|2007h . 
and in the MIR there are signs of recent or ongoing star formation 
jLacv et al.ll2005l : lYan et al]|2007h . 

At higher redshift there are samples of sub -mm selected 
QSO s with very high SFR (1000-3000 Mpyr ' ; iBertoldi et al] 
I2OO3I : iBeelen et alj l2006l : i Lutz etalJ |2007|, |2008|), but most of 
the high-z QSOs are u ndetected in the sub-mm range (~70 per 
cent; lOmont et alj|2003h . and the observed trend indicates that the 
starburst-to -AGN luminosity ratio decreases with i ncreasing AGN 
luminosity (H aas et al. 1 12003'; 'Maiol ino et idl bOOT*). thus suggest- 
ing that the observed correlation betwen star formation and AGN 

luminosity saturates at higher luminosities. 

The same pattern seems to apply to QS02s: ISturm et"ail 

( I2OO 6) find that IRS spectra of X -ray selected obscure d AGN show 
very weak or no PAH bands, and lPolletta eTal] ( |2008|) estimate that 
the star formation contributes <20 per cent to the bolometric lumi- 
nosity of a sample of obscured AGN with Lj^i^ > 10'" Lq. 

The ELAIS-IRS AGN comprise a very heterogeneous popu- 
lation, since it covers a wide range of PAH equivalent width and 
SFR. In some of the sources, the estimated starburst-induced IR 
luminosity (5FLir) calculated from the PAH luminosity is higher 
than the total IR luminosity (Ljr) obtained from the SED fitting. 
Since it is unlikely that Lpah/'S f Ar be lower for these sources (in- 
deed, PAH destruction by the AGN would lead to lower ratios) the 
most likely interpretation is that the SED of these AGN has a cold 
dust component which is not accounted for in the SED fitting - 
effectively underestimating Lir- because of lack of deep enough 
photometry in the FIR. Incidentally, MIPS70 and MIPS 160 upper 
limits for many ELAIS-IRS AGN cannot rule out a significant cold 
dust component in the FIR (see fits in Fig.[8]l. 

The PAH features are comparatively very weak in most 
ELAIS-IRS AGN, and in many of them (35 sources, mostly unob- 
scured AGN) only a loose upper limit of the (PAH-derived) 5 FLjr 
is available. To get a tighter estimate of its typical value, we mea- 
sure the strength of the PAH features in the averaged spectra of 
QSO and AGN2 sources. Assuming that the equivalent width of a 
PAH feature in the averaged spectrum (iiVVa) is similar to those of 
the individual sources, we can roughly estimate the typical PAH lu- 
minosity of the sources by multiplying EW^ times the luminosity of 
the continuum underlying the PAH feature. Since individual spectra 
are normalized at 7 \im before averaging them, the mean PAH lumi- 
nosity may be overestimated because the lower luminosity sources 
tend to show higher PAH equivalent widths. To minimize this bias, 
we split the QSO and AGN2 populations in two subsamples of high 
and low luminosity sources, based on their restframe 5.5-|,im lumi- 
nosity. We place the limit between populations at L55 = 4 x 10" Lq 
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Figure 22. Composite spectra of ELAIS-IRS AGN split into four groups 
based on their optical classification and 5.5-|.im restframe continuum lumi- 
nosity, (a) optical QSOs with vLy > 4 X lO" Lq. (b) optical QSOs with 
vLy < 4 X 10" Lq. (c) optical galaxies with vL,, > 4 X lO" Lq. (d) optical 
galaxies with vLy < 4 X lO" Lq. The dashed lines mark the position of the 
6.2, 7.7 and 1 1 .3 iim PAH features. 



to include the same number of sources in both subsamples. In the 
lower luminosity population the 6.2, 7.7 and 11.3 \im PAH features 
are detected in both QSO and AGN2 averages, while in the higher 
luminosity ones they all but disappear (Fig.l22t. 

TablellSlshows the estimated average PAH luminosity and as- 
sociated S FLin and SFR for the averaged spectra of the ELAIS-IRS 
AGN. In the lower luminosity QSOs the average PAH luminosity 
corresponds to an IR luminosity in the LIRG range, and accounts 
for 20-25 per cent of Ljr a s estimated from SED fitting, similar 
to the mean value found by IShi et al.l ( |2007|) in a sample of PG, 
2MASS and 3CR QSOs at lower redshift. The obscured AGN in 
the same luminosity range have PAH luminosities almost 3 times 
higher at 6.2 and 7.7[xm (which suggests star formation is respon- 
sible for about half of Lir in these sources), but they are similar 
to that of QSOs at 11.3 i^im. Estimated SFR is -50 Moyr"' for the 
lower-luminosity QSOs and ~100 Mgyr"' for the obscured AGN. 

In the higher-luminosity AGN the PAH features are too weak 
to meassure them even in the composite spectrum. Their equivalent 
width is at least 3 times lower in the high luminosity QSO compos- 
ite than in the lower luminosity one, and the resulting upper limits 
in LpAH constraint the mean SFR to < 100 Moyr"'. The implied 
IR luminosity from star formation accounts for <10 per cent of the 
estimated Lir. The high-luminosity obscured- AGN composite con- 
tains only five sources and has lower S/N than the others. Upper 
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limits for the 6.2 and 7.7 [xm features suggest SFR < 300 Mgyr"', 
and a starburst contribution to Lir of 10-15 per cent at most. 

To better constrain the amount of FIR emission in the QSOs, 
we have performed stacking analysis of the 70 and 160 |,im photom- 
etry. Fluxes for individual sources are estimated from the counts in 
the pixel containing the coordinates of the source, and their uncer- 
tainties are obtained from a fit to the histogram of (finite-valued) 
pixel values in the image. A mean sky level is computed for the im- 
age and subtracted from the source flux in order to ensure that the 
mean value for random positions in the image is zero. The mean 
flux X and dispersion ctq for the population are estimated from the 
individual measurements (x,, cr,) as the maximum-likelihood solu- 
tion for the distribution: 



p(x, cr) : 



1 



zexp\ 



2(0-2 + 0-2) j 



(8) 



A deta iled description of the whol e stack ing procedure can 
be found in ISerieant & Hatziminaogloul (l2009h . The mean fluxes 
obtained at 70 and 160 |xm are shown in Tablell4l If we assume the 
FIR SED is dominated by a starburst component with an M82-like 
SED, the average IR luminosity needed to reproduce the observed 
mean fluxes is ~10'2 Lq for the z < 1.5 sources, and ~5xl0'2 L© 
for those at z > 1.5. Similar results are obtained if we assume an 
Arp220 SED instead of M82 for the starburst component. 

The S FLiR extrapolated from the observed PAH luminosity is 
a factor 3 lower in the z < 1 .5 QSOs (and a factor < 10 lower in the 
z > 1.5 QSOs) than that obtained from the stacking analysis of 70 
and 160 \im photometry, suggesting that the average R\H-to-FIR 
luminosity ratio is significantly lower than in starburst galaxies. A 
likely interpretation is reduced PAH emission in star-forming re- 
gions in the vicinity of an AGN due to destruction of the PAH 
carriers by X-ray and far-UV photons from the AGN ( IVoitlll992l : 
iLe Floc'h et"ai]|200lh ISiebenmorgen et al.ll2004l) . 

Incidentally, the factor >3 reduction in EW-^ and in the PAH- 
to-FIR luminosity ratio between low- and high-luminosity QSOs 
suggests higher depletion rates of the PAH carriers in the most lu- 
minous AGN. Alternatively, it could indicate a saturation of the 
correlation between the luminosities of the starburst a nd AG N, as 
found in other samples of z < 0.5 QSOs (Shi et al] [jZOOTh and 
highe r-luminosity, high-z QSOs jMaiolino et al J 1 20071 ; Lutz et af] 
|2008[) . but results from the stacking analysis suggest the FIR lumi- 
nosity roughly scales with AGN luminosity. 

Estimated average luminosity of the 6.2 and 7.7 |.im PAH fea- 
tures in the low-luminosity composites are about 2-3 times higher 
in the obscured AGN than in the QSOs. This could indicate higher 
SFR in the obscured AGN, or, alternatively, higher abundance of 
PAH carriers due to absorption by the obscuring material of the en- 
ergetic photons from the AGN. However, the statistical significance 
of this result is not strong due to the small number of sources con- 
tributing to the obscured AGN average, and requires validation in 
larger samples. 

A selection bias favouring obscured sources with strong PAH 
emission could also be responsible for the increased PAH lumi- 
nosity, but we have not found any likely bias able to increase PAH 
luminosity in obscured AGN but not in QSOs. In example, the PAH 
features increase the observed flux in the selection band at z ~ 1, 
where we find roughly half of the obscured AGN, but these sources 
are bright enough to be selected even without contribution from 
PAH emission to the 15-|,im flux. Since the PAH features help to 
find a redshift estimate in many sources, this could play a role in 
the observed asymmetry, because most redshifts for AGN2 sources 



rely on the IRS spectrum, while most QSOs have optical spectro- 
scopic redshifts. But in AGN2 sources a reliable redshift estimate 
can usually be obtained from the (absorption) silicate feature alone, 
and indeed all sources with unreliable redshifts are classified as 
AGNl. In addition, a lower accuracy in redshifts obtained from the 
silicate feature relative to those from optical spectroscopy imply 
that PAH features may be blurred in the AGN2 composite, leading 
to an underestimation of their real strength. 

Another possible interpretation could be lower MIR luminos- 
ity in obscured AGN than in QSOs of the same intrinsic luminos- 
ity, as is expected if the putative torus is optically thick in the MIR; 
but this is not supported by measurements of the luminosity ra- 
tio between X-rays and MI R (Lutz e t al.ll2004l ; ISturm et al.ll2006l ; 
IPoIIetta et alJl2007l ; lHorst et aL200a) . 

6.4 Implications for the dust distribution 

A strong difference in SFR between optical type I and type 2 AGN 
has be en observed in a sample of QSOs and QS02s from the SDSS 
jKim et al. 2006). Using estimates from the [O n] 3727 A line, they 
find a SFR an order of magnitude higher in type 2 sources, and this 
is probably a lower limit, due to higher extinction in the optical. 

Such a large difference between type 1 and 2 sources can 
hardly be explained within the AGN unification scheme, and might 
indicate instead that the AGN classification into type I or type 2 
does not depend exclusively on the orientation of an AGN torus. 
Recent papers suggest that a fra ction of the obscuration in type 2 
AGN occurs in the host galaxy: iLacv et aL I ( l2007h find a correla- 
tion between silicate extinction and orientation of the plane of the 
galaxy, while PoUetta et al. (2008) suggest that foreground extinc- 
tion by cold dust is required to reproduce both the silicate feature 
and NIR continuum for almost half of the sources in a sample of 
MIR selected obscured AGN. 

If we assume that the extinction law for dust in the torus is 
similar to that found in local ULIRGs and in the Galactic Cen- 
tre, it reaches a minimum in opacity at ~7 |.im and increases to- 
wards shorter wavelengths; thus the optical depth is higher in the 
NIR than at 7 |,im. If the silicate absorption feature originates in the 
torus, a decrease in the viewing angle of the torus should produce 
a deeper silicate feature and weak er and steeper NIR continuum 
peaking at longer wave lengths (e.g. lEfstathiou & Rowan-RobinsonI 
1 19951 ; Ipritz et alj|2006l) . By contrast, if the silicate feature origi- 
nates in the host galaxy, a deep silicate profile could be observed in 
conjunction with an intense and flat NIR spectrum jPoUetta et al.l 

[2008;) . 

Fig.|16|compares the average spectrum of ELAIS-IRS AGN2 
sources with the average of sourc es that fit the torus model (T) 
or torus+cold absorber (T+C) from lPolletta et alj j2008l) . The NIR 
slope and strength of the silicate feature are better reproduced by 
the T model. 

Figure |23] shows the averaged spectrum of the ELAIS-IRS 
AGN separated in four groups attending to their measure of 5sii. 
In sources with 5sii > (mostly QSOs) the NIR slope is small and 
independent of 5sii, but in sources with 5sii < the NIR slope is 
substantially higher, particularly for the subsample with strongest 
silicate absorption. This confirms the existence of a correlation be- 
tween NIR slope and silicate absorption, as expected if the obscu- 
ration originates mainly in the AGN torus. 

To search for sources that depart from this trend, we calculated 
the spectral index a of the continuum between 2.2 and 7 |xm rest- 
frame for individual spectra (Fig. 1241 ). Most QSOs yield a < 1.5, 
while most AGN2 show a > 1.5. There is a clear trend, even if the 
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Figure 23. Composite spectra of ELAIS-IRS AGN separated into 4 groups 
based on tlieir measure of silicate strength: S^u < -0.5 (solid line), -0.5 < Ssii 
< (dotted line), < S^n < 0.15 (dashed line), and S^ii > 0.15 (dot-dashed 
line). 
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Figure 24. Relation between silicate strength and spectral index of the con- 
tinuum between 2.2 and 7 [im restframe for ELAIS-IRS sources. Symbols 
represent: unobscured AGN (plus signs), obscured AGN (diamonds), and 
starbursts (asterisks). 



dispersion is high, towards a simultaneous increase of a and 5,0 in 
AGN, but not in starbursts. Note that EIRS-21, which is suspected 
of harboring an obscured AGN due to its strong MIR continuum, 
departs strongly from the AGN population. This suggests that the 
observed extinction is better reproduced by a screen of cold dust 
instead of a thick AGN torus. Incidentally, it is undetected down 
to 1 X 10""' erg cm"' s"' in [2-10] keV band images from XMM 
(I. Valtchanov priv. comm.), so the presence of a dust enshrouded 
AGN is questionable. 

The only AGN2 source that separates slightly from the general 
trend is EIRS-52, but its spectrum is still better described by the 
torus model; therefore, a dusty torus is favored as the source of the 
obscuration in the ELAIS-IRS AGN. 

The lack of T-l-C sources in ELAIS-IRS is probably caused by 
the selection criteri a, which prevent se lection of highly obscured 
AGN as those from lPolIetta et al.l ( |2008|) . 



7 CONCLUSIONS 

We have analyzed the IR properties of a sample of 70 MIR se- 
lected LIRGs and ULIRGs at 0.4 < z < 3.2. The MIR selection 
favours starburst-dominated sources in redshift intervals in which 
the strongest PAH bands overlap with the passband of the selec- 
tion filter, namely, z ~ 1 at 15 |xm and z ~ 1.8 at 24 |xm. For the 
AGN-dominated sources, the redshift distribution is wider because 
their IR SED has smoother features. The relatively shallow depth 
of the optical data in the ELAIS-IRS sample and the requirement 
of optical photometric or spectroscopic redshifts for source can- 
didates favours the selection of weakly obscured sources, contrar- 
ily to other MIR high-z ULIRG samples where highly obscured 
sources are preferred, and thus turns the ELAIS-IRS sample into 
a valuable complement for the study of the high-z ULIRG popula- 
tion. 

Redshift estimates for the ELAIS-IRS sources from their MIR 
spectrum are reliable for obscured or star-forming sources, and un- 
reliable in many unobscured AGN. This behaviour is the opposite 
to that obtained through optical spectroscopy, and thus guarantees 
that reliable redshift estimates can be obtained by one or other 
means for most IR sources. 

A combination of MIR diagnostic criteria efficiently separates 
starbursts, obscured and unobscured AGN, with a good correlation 
to the optical classification. The selection criteria of the ELAIS- 
IRS sample favour the selection of starburst-dominated sources at 
0.6 < z < 1.2, yet most galaxies selected in this range are AGN- 
dominated. This suggests that the population of 5i5 > 1 mJy sources 
at z ^ 1 is dominated by AGN. 

The averaged spectrum of ELAIS-IRS starbursts is similar to 
that of local ULIRGs, but the observed depth of the silicate fea- 
ture, strength of PAH features and PAH flux ratios are more alike to 
those found in lower-luminosity starburst galaxies. This behaviour 
has already been observed in higher-redshift starburst samples and 
could be related to star formation extending over larger regions or 
a higher gas-to-dust ratio in the galaxy. 

The correlation between MIR and FIR SEDs allows the es- 
timation of IR luminosities from the MIR spectrum with a factor 
3 uncertainty. If MIPS70 and/or MIPS 160 photometry is available 
the uncertainty reduces further and becomes dominated by uncer- 
tainties in the photometry and in the SED models. 

The mean SFR in ELAIS-IRS QSOs estimated from PAH lu- 
minosity is 50-100 Moyr"', but its implied FIR emission is 3-10 
times lower than that obtained from stacking analysis in the 70 and 
160 \im bands, thus suggesting the PAH emission in substantially 
depleted in most QSOs. 

A lower PAH-to-continuum luminosity ratio is found in the 
most luminous QSOs relative to lower-luminosity ones, suggesting 
a saturation of the correlation between star formation and AGN ac- 
tivity found in other samples of lower-luminosity sources or, alter- 
natively, higher depletion rates of the PAH carriers in the galaxies 
hosting the most powerful AGN. 

In the obscured AGN, the mean SFR is 2-3 times higher than 
in QSOs of the same luminosity. If confirmed in larger samples, 
this could indicate a connection between dust obscuration and star 
formation in the host galaxy. However, the observed correlation be- 
tween silicate absorption and the slope of the NIR-to-MIR contin- 
uum is consistent with the hypothesis of the obscuring dust being 
located in a dusty torus. 
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Table 1. IRS observations log 



ID 




ELAIS name 


R.A." 


Dec." 






PA. 


0b.s 


date 




Exposure time (s) 










[h.m..s.] 


[d.m.s] 


[mJy] 




[deg] 








SL2 


SLl 


LL2 


LLl 


EIRS 


1 


ELAISC 15J161401.0+544733 


16hl4m30.9s 


+54d52ml5s 


0.79 + 


10 


147.1 


2005 


01 


08 


967 


967 


1950 


1950 


EIRS 


2 


ELAISC15J161350.0+542631 


16hl4ml9.3s 


+54(13 1ml Is 


0.80 + 


16 


146.9 


2005 


01 


08 


1451 


967 


1950 


1950 


EIRS 


3 


ELAISC15J163636.2+411049 


16h36m42.9s 


+41dl7m04s 


0.83 + 


14 


116.2 


2005 


02 


14 


967 


967 


1950 


1950 


EIRS 


4 


ELAISC 1 5 J003640-433925 


00h36mll.8s 


-43d35m46s 


0.98 ± 


17 


50.0 


2005 


06 


08 


967 


967 


1950 


1950 


EIRS 


5 


ELAISC15J003408-4310I1 


00h34m4l.4s 


-43dl2m3ls 


0.99 ± 


13 


216.2 


2004 


11 


16 


967 


967 


1950 


1950 


EIRS 


6 


ELAISC15J163721. 3+41 1503 


16h37m29.2s 


+41d21ml3s 


1.00 + 


14 


118.4 


2005 


02 


12 


484 


484 


975 


975 


EIRS 


7 


ELAISC15J163655.8+405909 


16h37m03.2s 


+41d05m23s 


1.02 + 


14 


117.4 


2005 


02 


13 


484 


484 


975 


975 


EIRS 


8 


ELAISC 15J163455.0+4 12211 


16h35m04.3s 


+41d28ml6s 


1.04 + 


15 


120.9 


2005 


02 


08 


484 


484 


975 


975 


EIRS 


9 


ELAISC 1 5 J 160733.7+534749 


16h08mOI.6s 


+53d52m41s 


1.04 + 


18 


144.5 


2005 


01 


09 


484 


484 


975 


975 


EIRS 


10 


ELAISC15J16151 1.2+550627 


16hl5m41.3s 


+55dllm06s 


1.06 + 


12 


147.5 


2005 


01 


08 


484 


484 


975 


975 


EIRS 


11 


ELAISC 1 5 J 163739.3+414348 


16h37m47.5s 


+41d49m59s 


1 .06 ± 


11 


118.4 


2005 


02 


12 


484 


484 


975 


975 


EIRS 


12 


ELAISC15J161229.0+542832 


16hl2m57.7s 


+54d33ml9s 


1.06 + 


17 


145.8 


2005 


01 


09 


484 


484 


975 


975 


EIRS 


13 


ELAISC15J163954.5+411109 


16h40m02.3s 


+41dl7m20s 


1.07+0 


16 


118.1 


2005 


02 


13 


484 


484 


975 


975 


EIRS 


14 


ELAISC15J163536.6+404754 


16h35m44.4s 


+40d54m05s 


1.12 + 


12 


118.1 


2005 


02 


12 


484 


484 


975 


975 


EIRS 


15 


ELAISCI5J161332.4+544830 


16hI4m02.5s 


+54d53ml3s 


1.12 + 


19 


146.9 


2005 


01 


08 


484 


484 


975 


975 



"coordinates of the optical ID of the 15-|im ELAIS source jGonzalez-Solares et alj200^ . 



Table 2. Optical and Near infrared photometry 



ID 


U 


S 




f 


Z 


J 


H Ks 


EIRS-1 


19.39 ± 0.02 


20.28 ± 0.03 


19.64 ± 0.04 


1941 +0.11 


18.95 ± 0.02 


18.17 + 0.02 


16.42 ±0.01 


EIRS-2 


20.51 ± 0.04 


23.08 ± 0.05 


22.14 + 0.06 


21.25+0.12 


20.40 ± 0.04 


19.10 + 0.07 


17.11 ±0.03 


EIRS-3 




23.22 ± 0.10 


22.27 ± 0.05 


21.01 ±0.10 


20.40 ± 0.05 






EIRS-4 
















EIRS-5 
















EIRS-6 


20.51 ±0.03 


20.77 ± 0.02 


20.59 ± 0.04 


20.14 + 0.06 


19.97 + 0.04 






EIRS-7 


22.63 +0.10 


22.87 ± 0.04 


22.92 ± 0.08 


23.02 ± 0.23 








EIRS-8 


21.94 + 0.05 


22.39 ± 0.03 


21.80 + 0.06 


20.90 + 0.10 


20.20 ± 0.05 






EIRS-9 






22.71 ± 0.07 


21.88 + 0.06 


21.11 + 0.12 






EIRS- 10 


19.07+0.02 


19.58 ± 0.02 


19.07 ± 0.02 


19.00 ± 0.02 


18.72 + 0.02 






EIRS-1 1 


18.22 + 0.02 


19.06 + 0.02 


18.79 ± 0.02 


18.51 ±0.02 


18.55 ± 0.03 






EIRS-12 


20.63 ± 0.05 


21.47+0.02 


21.14 + 0.02 


20.79 + 0.11 


20.25 ± 0.04 


20.23 ±0.13 


19.46 ±0.15 


EIRS- 13 




23.57 ± 0.07 


22.76 + 0.06 


22.10±0.11 


20.49 ± 0.06 






ElRS-14 






23.46 + 0.14 


21.78±0.12 


21.18 + 0.11 






EIRS-15 


22.49 + 0.10 


22.72 + 0.04 


21.74 + 0.05 


20.49 ±0.11 


19.68 + 0.02 


18.08 ±0.02 


16.41 ±0.01 



This paper has been typeset from a TgX/ LSTgX file prepared by the author. 



Table 3. mid- and far-IR photometry 



ID 




S4.5 


^5.8 




Sl5 


^24 


-^70 


-^'160 




li'Jy] 


liJyl 


[liJy] 


[liJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


ElRS-1 


297.0 ± 1.9 


366.0 + 2.4 


471 +6 


789 + 5 


0.79 + 0.10 


2.48 + 0.02 






ElRS-2 


128.0+ 1.6 


113.0+ 1.4 


64 + 5 


84 + 4 


0.80 ±0.16 


0.94 + 0.02 


164 + 1.0 


70.2 + 5.1 


EIRS-3 


63.4 + 1.1 


40.7 + 1.2 






0.83 ±0.14 








EIRS-4 


347.0 + 2.0 


524.0 ± 2.6 


745 + 7 


1090 ± 7 


0.98 ±0.17 


4.33 + 0.04 


33.1 + 1.7 




ElRS-5 


276.0 + 2.2 


359.0 ± 2.7 


458 + 8 


707 ± 8 


0.99 ±0.13 


2.18+0.03 






ElRS-6 


31.2 + 0.9 


41.6 ± 1.2 


56 + 4 


140 ±4 


1.00 ±0.14 


0.48 + 0.02 






ElRS-7 


27.6 + 0.5 


52.6 ± 1.0 


121 + 3 


399 ±4 


1.02 ±0.14 


2.39 ± 0.02 


13.0+1.2 




ElRS-8 


205.0+ 1.8 


247.0 ± 1.6 


345 +6 


516 + 4 


1.04 ±0.15 


1.40 ± 0.02 






ElRS-9 


98.6 + 1.3 


77.7 ± 1.5 


96 + 5 


163 + 5 


1.04 ±0.18 


1.91 ± 0.02 


36.1 + 1.3 




EIRS-10 


560.0 + 2.9 


722.0 ± 3.6 


823 + 7 


1100 + 6 


1.06 ±0.12 


2.73 ± 0.02 






EIRS-1 1 


171.0+ 14 


267.0 ± 1.7 


412 + 5 


693 +4 


1.06 ±0.11 


2.34 + 0.02 






EIRS-12 


27.1 +0.7 


45.0 ± 1.1 


102 + 4 


248 + 4 


1.06 ±0.17 


2.11 +0.02 






EIRS-13 


99.4 + 1.2 


81.2 ± 1.2 


68+4 




1.07 ±0.16 


0.48 + 0.02 






EIRS- 14 


107.0+ 1.0 


82.8 + 1.0 


109 + 3 


133 + 3 


1.12±0.12 


1.30 + 0.02 


22.7 + 14 


66.1 +5.5 


EIRS-15 


264.0+ 1.8 


230.0+ 1.7 


336 + 6 


347 + 4 


1.12 ±0.19 


141 +0.02 




47.3+5.9 



Infrared-luminous galaxies at z 0.5 



Table 4. Library of mid-IR spectra 



ID 




ID ; 


ID 








ULIRGs; Armus ct al. 2007 






IRAS 05189-2524 
IRAS 14348- 1445 
Arp220 
UGC5I0I 


0.043 
0.083 
0.018 
0.039 


IRAS 08572-1-3915 0.058 
IRAS 15250+3609 0.055 
Mrk231 0.042 

Starbursts; Brandl et al. 2006 


IRAS 12112+0305 
IRAS 22491-1808 
Mrk273 


0.073 
0.078 
0.038 


Mrk52 
NGC 3628 
NGC7714 


0.007 
0.003 
0.009 


NGC 1222 0.008 
NGC 520 0.008 


NGC 2146 
NGC 7252 


0.003 
0.016 



Table 5. Redshifts 



ID 






opt. cla.ss'' 


;lypc 




Reference 


EIRS 


-1 


0.387 


QSO 


zspec 


A 


Perez-Foumon (priv. comm.) 


EIRS- 


-2 


1.154 


galaxy 


zIRS 


A 




EIRS- 


-3 


0.676 


galaxy 


zIRS 


A 




EIRS 


-4 


1.181 


QSO 


zspec 


A 


La Franca et al. (2004) 


EIRS 


-5 


1.065 


QSO 


zspec 


A 


La Franca et al. (2004) 


EIRS 


6 


2.356 


QSO 


zspec 


A 


Perez-Fournon (priv. comm.) 


EIRS 


-7 


2.592 


galaxy 


zspec 


A 


Swinbank et al. (2004) 


EIRS 


-8 


0.884 


galaxy 


zIRS 


B 




EIRS 


■9 


0.609 


galaxy 


zIRS 


A 




EIRS- 


-10 


1.265 


QSO 


zIRS 


B 




EIRS 


11 


1414 


QSO 


zspec 


A 


SDSS 


EIRS 


-12 


2.024 


galaxy 


zIRS 


A 




EIRS 


-13 


1.091 


galaxy 


zIRS 


A 




EIRS 


-14 


0.619 


galaxy 


zIRS 


A 




EIRS 


-15 


0.827 


galaxy 


zIRS 


A 





'^Classification of the source based on its optical SED or spectrum 
''Reliability of the redshift estimation. A = reliable, B = unreliable 



Table 6. Monochromatic fluxes and luminosities 



ID 




/,. 2.2 fim 


fy 5.5 fim 


fr 7 nin 


fy 10 nm 


fy 15 ^m 


vLy 2.2 nm 


vLy 5.5 fim 


vLy 1 pm 


vLy 10 pm 


vLy 15 pm 






[mJyl 


[mJyl 


[mJy] 


[mjyj 


[mJy] 


[10'»Lo] 


llO'°Lol 


[10'°Lo] 


[10'»Lol 


[10'°Lo] 


EIRS 


-1 




0.67+0.10 


1.0 + 0.1 


14 + 0.1 


24 + 0.1 




3 + 


4 + 


3 + 


4 + 


EIRS 


-2 


0.104 + 0.004 


0.18+0.07 


0.5+0.1 


0.5+0.2 


1.0 + 0.2 


13.0 + 0.5 


9 + 3 


21 + 3 


13 + 4 


19 + 3 


EIRS 


-3 


0.061 + 0.002 


<0.20 


<0.2 


<0.2 


<0.2 


2.6 + 0.1 


<3 


< 2 


< 1 


< 1 


EIRS 


-4 


0.575 ± 0.008 


1.56 + 0.10 


2.0 + 0.1 


3.9 + 0.2 


6.3 ± 0.2 


75.0+ 1.0 


81+5 


82 + 4 


111 +5 


120 + 3 


EIRS 


-5 


0.362 + 0.006 


1.08+0.09 


1.2 + 0.1 


1.6 + 0.2 


3.0 + 0.2 


38.5+0.6 


45+3 


41 + 3 


37 + 4 


47 + 2 


EIRS 


-6 


0.117 + 0.009 


<0.40 


0.3 + 0.2 


0.9 + 04 


1.0 + 0.2 


55.1 +4.2 


<75 


50 + 24 


94 + 38 


69 + 10 


EIRS 


-7 


0.387 ± 0.009 


2.33 + 0.25 


2.7+0.2 




8.3 + 1.3 


215.2 + 5.0 


518 + 54 


466 + 26 




680 + 102 


EIRS 


-8 


0.230 + 0.003 


0.66 + 0.10 


0.9 + 0.1 


0.9 + 0.2 


1.8 + 0.2 


16.9 + 0.2 


19 + 2 


21+2 


14 + 3 


19 + 1 


EIRS 


-9 


0.099 ± 0.001 


0.35+0.10 


0.9 + 0.1 


04 + 0.1 


1.8 + 0.1 


34 ± 0.0 


4 + 1 


9 + 


3 + 1 


8 + 


EIRS 


-10 


0.760 ± 0.010 


1.37+0.11 


1.9 + 0.1 


2.8 + 0.2 


4.4 + 0.4 


1134+ 1.5 


81+6 


89 + 5 


90 + 6 


96 + 7 


EIRS 


11 


0.357 + 0.009 


1.12+0.17 


1.5+0.1 


2.9 + 0.2 


4.7 ± 0.7 


66.0+ 1.7 


82 + 12 


88 + 8 


119 + 6 


126 + 19 


EIRS 


-12 


0.159 + 0.009 


1.40 + 0.16 


1.9 + 0.2 


1.8 + 0.2 


4.1 +0.6 


57.2 + 3.2 


202 + 22 


214 + 23 


143 + 14 


216 + 32 


EIRS 


-13 


0.080 + 0.003 


<0.21 


0.3+0.1 


<0.5 


<0.4 


8.9 + 0.3 


<9 


10 + 4 


< 11 


<6 


EIRS 


-14 


0.107 + 0.001 


<0.29 


0.7+0.1 


0.3 + 0.1 


1.1 ±0.1 


3.8 + 0.0 


<4 


7 + 1 


2+ 1 


5+0 


EIRS 


-15 


0.248 ± 0.004 


0.51 +0.10 


0.6 + 0.1 


0.8 + 0.2 


1.6 + 0.1 


15.9 + 0.3 


12 + 2 


12 + 2 


11+2 


14+ 1 



A. Herndn-Caballero et al. 



Table 7. PAH fluxes and luminosities. Silicate absorption 



ID l0'5x/62 EWf,2 Lf,2 l0i5x/77 EW^^ L^^ T9.7 

[ergcm-2 s"'] [|im] [lO'^Lo] [erg cm-^ s"' ] [|im] [lO'" Lo] 



EIRS 


1 


<4.6 


<0.11 


<0.06 


14.0 ±2.0 


0.38 ± 0.06 


0.19 + 0.03 


0.00 + 0.11 


0.05 


EIRS 


2 


9.6 ± 0.9 


1.43 ±0.14 


1.89 ±0.18 


40.0 ± 1.5 


5.70 ± 0.21 


7.88 ± 0.29 


0.61 ±0.61 


-0.58 


EIRS 


3 


6.7 ± 2.2 


5.86 ± 2.00 


0.35 ± 0.12 


14.3 ± 2.6 


13.61 ±2.60 


0.74 ± 0.14 


-1.82 ± 2.09 


1.93 


EIRS 


4 


<2.5 


<0.04 


<0.52 


<2.5 


<0.05 


<0.52 


-0.23 ± 0.08 


0.17 


EIRS 


5 


5.1 ± 1.5 


0.12 ± 0.04 


0.81 ±0.24 


5.1 ± 1.5 


0.16 ± 0.05 


0.83 ± 0.25 


0.08 ±0.15 


-0.07 


EIRS 


6 


<2.8 


<0.31 


<3.20 


<2.4 


<0.34 


<2.76 


-0.22 ± 0.50 


0.28 


EIRS 


7 


< 2.7 


< 0.05 


< 3.90 


< 2.2 


<0.05 


<3.20 


0.95 + 0.25 


-0.76 


EIRS 


8 


<4.4 


< 0.15 


< 0.45 


<6.9 


< 0.27 


<0.70 


0.48 ± 0.25 


-0.51 


EIRS 


9 


22.0 ± 1.7 


2.00 + 0.15 


0.89 ± 0.07 


82.9 ± 2.7 


7.32 + 0.24 


3.34 + 0.11 


0.77 ± 0.38 


-0.58 


EIRS 


10 


< 3.5 


< 0.07 


<0.86 


< 3.0 


< 0.07 


<0.73 


-0.05 + 0.11 


0.13 


EIRS 


11 


7.8 ± 2.6 


0.20 ± 0.07 


2.56 ± 0.87 


6.7 ± 2.3 


0.19 ± 0.07 


2.20 ± 0.76 


-0.31 ±0.19 


0.29 


EIRS 


12 


5.9 ± 2.4 


0.17 ± 0.07 


4.74 ± 1.88 


13.4 ±3.1 


0.42 ± 0.10 


10.70 ± 2.48 


0.95 ± 0.21 


-0.77 


EIRS 


13 


5.3 ± 1.6 


1.21 ± 0.38 


0.90 ± 0.28 


21.2 ±2.4 


5.49 ± 0.63 


3.64 ± 0.42 


0.80 ± 1.98 


-00 


EIRS 


14 


25.0 ± 2.9 


3.11 ±0.37 


1.05 ±0.12 


64.5 ± 3.5 


8.47 ± 0.47 


2.71 ±0.15 


0.84 ± 0.60 


-0.38 


EIRS 


15 


2.3 ± 1.1 


0.10 ± 0.05 


0.20 ± 0.09 


10.6 ± 2.2 


0.59 ± 0.12 


0.91 ± 0.19 


0.13 ± 0.28 


-0.02 



Table 8. Diagnostic criteria 



Criterion 


Expression 


AGNl 


AGN2 


SB 


CI 


£W,,2 


< 0.2 


<0.2 


>0.2 


C2 




>0.1 


>0.1 


<0.1 


C3 


'TDR 


<15% 


<15% 


> 15% 


C4 


Ucn 


<0.15 


0.15-1.0 


> 1.0 


C5 


5sil 


>0 


<0 


<0 


C6 


L15/L10 


<2 


2-4 


>4 



Table 9. IR vs optical classification for the Library 



Optical class 


N AGNl 


N AGN2 


NSB 


Syl 


31 


2 





Sy 1/HlI 


3 


3 


1 


Sy2 


2 




1 


Sy2/HII 


1 


1 


15 


KI 








7 


ULIRG/LINER 





3 


12 


ULIRG/HII 








12 



Table 10. Diagnostic results for ELAIS-IRS 



ID 


z 


EW62 [nm] 


CI 




C2 


rpDR m 


C3 


Hen 


C4 




C5 




C6 


Cflnal 


EIRS-1 


0.387 


< 0.112 


1.2 


0.486 


1,2 


3.9 


1.2 


0.260 


2 


0.05 


1 


1.777 


1 


1 


EIRS-2 


1.154 


1.434 


S 


0.030 


s 


57.1 


S 


8.010 


S 


-0.58 


2,S 


2.068 


2 


S 


EIRS-3 


0.676 


5.861 


S 


< 0.268 


S 


57.5 


S 


13.880 


S 


1.28 


1 


0.721 


1 


S 


EIRS-4 


1.181 


< 0.040 


1,2 


0.225 


1,2 


0.8 


1,2 


< 0.046 


1 


0.16 


1 


1.614 


1 


1 


EIRS-5 


1.065 


0.125 


1,2 


1.024 


1,2 


3.7 


1,2 


0.310 


2 


-0.07 


2,S 


1.875 


1 


2 


EIRS-6 


2.356 


< 0.306 


1,2 


< 0.558 


S 


0.0 


1,2 


< 0.340 


1 


-0.26 


2,S 


1.101 


1 


1 


EIRS-7 


2.592 


< 0.052 


1,2 


0.316 


1,2 


0.0 


1,2 


< 0.052 


1 


-0.76 


2,S 






2* 


EIRS-8 


0.884 


< 0.152 


1,2 


0.528 


1,2 


1.1 


1.2 


< 0.274 


1 


-0.51 


2.S 


2.011 


2 


1 


EIRS-9 


0.609 


2.004 


S 


0.046 


S 


46.0 


S 


6.560 


S 


-0.58 


2,S 


4.036 


S 


S 


EIRS- 10 


1.265 


< 0.066 


1.2 


0.899 


1,2 


0.0 


1.2 


< 0.066 


1 


0.20 


1 


1.588 


1 


1 


EIRS-1 1 


1.414 


0.199 


1.2 


0.654 


1,2 


3.3 


1.2 


0.300 


2 


0.17 


1 


1.599 


1 


1 


EIRS-12 


2.024 


0.165 


1,2 


0.117 


1,2 


4.9 


1,2 


0.260 


2 


-0.77 


2,S 


2.264 


2 


2 


EIRS-13 


1.091 


1.208 


S 


< 0.084 


S 


51.5 


S 


5.920 


S 


-0.80 


2,S 


> 1.069 


S 


S 


EIRS-14 


0.619 


3.113 


s 


< 0.048 


S 


45.7 


S 


7.690 


s 


-0.38 


2,S 


3.486 


2 


S 


EIRS-15 


0.827 


0.103 


1,2 


0.878 


1,2 


11.6 


1,2 


0.900 


2 


-0.02 


2,S 


1.893 


1 


2 



'Assuming the unknown diagnostic from L15/L10 agrees with that of S^i\. 
■•"Source classifies as both AGNl and SB, but priority is given to the later. 



Infrared-luminous galaxies atz~ 0.5-3 3 1 



Table 11. IR vs optical classification for ELAIS-IRS 



Optical class 


N AGNl 


N AGN2 


NSB 


QSOs 


33 


5 


1 


Galaxy 


7 


12 


10 



Table 12. Recent measurements of LpWiiR from the Literature 



Reference 


Galaxy type 


PAH teature(s) 


Result 


Rigopoulou et al. (1999)' 


starburst 


7.7 (un 


L77/L1R = 0.0081 




ULIRG starburst 


7.7 ^im 


L77/I1R = 0.0055 


Liitz et al. (2003)'' 


starburst 


7.7 ^im 


LvT/ilR = 0033 


Spoon ct al. (2004)" 


normal and starburst 


6.2 pm 


L62/i-IR = 00034 


Smith et al. (2007)' 


normal and starburst 


6.2 pm 


L62/L,R= 0.011 




normal and starburst 


7.7 pm 


Ln/Lm = 0.041 


Farrah et al. (2008)' 


starburst 


6.2+11.3 pm 


L62+113/AR = 0.038 


This work 


ULIRG starburst 


6.2 pm 


L62/L1R = 0.012 




ULIRG starburst 


7.7 pm 


L77/ilR = 0.038 




ULIRG starburst 


11.3 [xm 


Lll3/ilR = 0.011 



^Estimates underlying continuum by interpolation in PAH-free intervals. 

'^Uses fv in the peak of the feature. Converted to luminosity by Shi et al. (2007) assuming a Drude 
profile with FWHM = 0.6 urn. 

'^Estimates continuum by simultaneous fitting of the spectiiim to a Lorentzians+continuum model. 



Table 13. IR Luminosity and SFR 



ID S70 S16O logiiR/Lo SED' PAH" SFLir/Li^' SFR" 

[mjy] [mjy] [Moyr"'] 



EIRS-1 






10.8+0.3 


QSO high 


6.2 11.3 




< 0.44 


< 5 


ElRS-2 


16.4 ± 1.0 


70.2 ± 5.1 


12.4 + 0.1 


Dale26 


6.2 7.7 11 


3 


0.69 ± 0.03 


332+ 12 


ElRS-3 






11.0 + 0.3 


Seyfert 2 


6.2 7.7 11 


3 


2.03 + 0.34 


40 + 6 


EIRS-4 


33.1 ± 1.7 




12.5+0.5 


Strong FIR QSO 


6.2 7.7 




< 0.08 


<47 


ElRS-5 






11.6 ±0.3 


QSO low 


6.2 77 11 


3 


1.47 ±0.19 


115 + 15 


EIRS-6 






12.1 ± 0.3 


QSO high 


6.2 7.7 




< 1.20 


<287 


EIRS-7 


13.0 ± 1.2 




13.0 ± 0.5 


Mrk231 


6.2 7.7 




<0.20 


<346 


EIRS-8 






11.3 ± 0.3 


QSO low 


6.2 7.711 


3 


< 1.25 


<43 


EIRS-9 


36.1 ± 1.3 




11.9 ±0.2 


NGC 6240 


6.2 77 11 


3 


0.88 ± 0.03 


124±4 


ElRS-10 






12.1 ±0.3 


Weak FIR QSO 


6.2 7.7 




< 0.31 


<77 


EIRS-1 1 






12.0 + 0.3 


QSO norm 


6.2 7.7 11 


3 


1.06 + 0.22 


223 + 45 


EIRS-12 






13.1 +0.3 


M82 


6.2 7.7 




0.25 ± 0.06 


576+ 143 


ElRS-13 






12.0 + 0.3 


GL12 


6.2 7.7 11 


3 


0.99 + 0.10 


173 + 17 


ElRS-14 


22.7 ± 1.4 


66.1 ±5.5 


11.9 + 0.1 


Dale26 


6.2 77 11 


3 


0.86 ± 0.05 


119±6 


EIRS-15 




47.3 ± 5.9 


11.3 ± 0.3 


Weak FIR QSO 


6.2 7.7 11 


3 


1.12 ±0.26 


42 ±9 



^template SED that best fits the MIR and FIR photometry. 
*'PAH feature(s) used for star formation rate estimation, 
'^ratio of star formation induced IR luminosity to total IR luminosity, 
•^star formation rate as estimated from PAH luminosity 



Table 14. 70 and 160fim stacking analysis of unobscured AGNs 



redshift 


filter 


mean flux 


L,R M82» 


LiR Arp220^ 






[mJy] 


[Lol 


[Lol 


z < 1.5 


70 pm 


11 + 2 


8.5x10" 


84x10" 


(22 sources) 


160 pm 


27 + 9 


1.3xl0'^ 


8.9xl0" 


z> 1.5 


70 pm 


5.3+ 1.1 


4.4x10'^ 


8.3xl0'2 


(18 sources) 


160 (xm 


28 ±7 


6.3x10'^ 


5.3x10'^ 



^Infrared luminosity of the starburst component required to reproduce ob- 
served FIR flux, assuming M82 or Arp220 SED. 



A. Herndn-Caballero et al. 



Table 15. SFR and S FLm in averaged spectra 



Group 


PAH band 




<^PAH > 


(^^■ilR> 


<SFR) 






[urn] 


[Le] 


[Le] 


[Moyr-1] 


QSOs 


6.2 ^rm 


<0.01 


<7.1 X 10' 


<5.9 X 10" 


<100 


high lum. 


7.7 |im 


<0.03 


<1.7 X 10'" 


<4.5 X 10" 


<80 


(18 sources) 


11.3 nm 


<0.01 


<3.5 X 10" 


<3.1 X 10" 


<53 


QSOs 


6.2 nm 


0.03 


3.3 X lO' 


2.7 X 10" 


46 


low kim. 


7.7 [Mir 


n.n'i 


s.i X in'' 


2,1 X Id" 




(18 sources) 


1 1 ..1 [inr 


0,(J.i 


.i.4 X 1(J'' 


.i.Ox 10" 




Obsc. AGN 


6.2 \im 


<0.04 


<2.5 X 10"' 


<2.1 X 10'^ 


<352 


high lum. 


7.7 ]im 


<0.10 


<5.6 X 10'" 


<1.5 X 10'2 


<257 


(5 sources) 


1 1.3 \im 


<0.20 


<3.8 X 10'" 


<3.5 X 10'2 


<589 


Obsc. AGN 


6.2 \im 


0.10 


8.1 X 10' 


6.7 X 10" 


114 


low lum. 


7.7 [tm 


0.34 


2.6 X 10'" 


6.8 X 10" 


121 


(6 sources) 


11.3 (un 


0.07 


3.6 X 10' 


3.3 X 10" 


54 



